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THEORETICAL CONSIDERATIONS ON POTENTIATION 
IN DRUG INTERACTION 


H. D. LANDAHL 
COMMITTEE ON MATHEMATICAL BIOLOGY 
THE UNIVERSITY OF CHICAGO 


To account for some of the more important aspects of drug interaction 
we shall consider a model which can also account for certain general 
properties of the action of a single drug. A simple model in which there 
may be enzymatic detoxification of a drug is studied theoretically. The 
relation between time for appearance of an effect due to the drug and the 
size of the dose is found to contain the same parameters as the relation 
between the effectiveness of paired doses and the interval of time be- 
tween doses. A similar situation holds when the drug is given at a con- 
stant rate. 

When two drugs are administered together, their effect will depend on 
the manner of interaction, how much of each drug is given, which is given 
first, and on the interval of time between each administration. A number 
of plausible types of interaction is considered theoretically in terms of 
the model, analytical expressions being given for a number of cases. 
The interaction may be synergistic or antagonistic. In the former case 
the potentiation may be more than or less than additive depending on the 
order of delivery and on the time between injections. Methods for the 
estimation of the parameters from data are discussed. 


The interaction between drugs is of considerable importance in 
chemotherapy, in the interaction of antibiotics and in the inter- 
action of toxic substances, such as are used for insecticides (Jawetz 
and Gunnison, 1953; Loewe, 1955; Frawley et al., 1957). Two 
drugs may act together synergistically or they may act antagonisti- 
cally. The former action we shall refer to as positive potentiation; 
the latter as negative potentiation. It is the purpose here to con- 
sider the problem in terms of a simple model in which temporal 
aspects are included. The model for time of effect is not as gen- 
eral in some respects as that considered by J. F. Yeager and S. C. 
Munson (1945). Similarly the models for interaction of drugs with 
enzymes is not as general as that considered by F. H. Johnson, 


i 
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H. Eyring, and W. Kearns (1948; cf. also Johnson, Eyring, and 
Polissar, 1954, chap 10). 


I. A simple model to represent the time of effect as a function 
of dose and the influence of injection rate 


Let a drug X produce an effect of magnitude z. We assume 2 
varies with time according to 
dz 
— =Ar-az, (1) 
dt 
where z is the concentration of X and where A and a are constants. 
If z represents the concentration of a substance then the term az 
measures the rate at which the substance is metabolized and ex- 
creted. If initially z is negligible, we obtain 


t 
oe Aew st [ d*’ x(t’) dt’, (2) 
0 


A being the proportionality constant. 
We assume that an enzyme E£ is responsible for the elimination 
of 2 so that we have the following reaction: 


Re Be ey 


where K is the ratio (X)(E)/(XE), and X, is a detoxified form of 
X. For simplicity we assume the reaction rate m to be much slower 
than the others. The above equation gives: 


(3) 


where E’, is the total enzyme concentration and 6 = mE,/K. There- 
fore 2 is given by 


x 
Kot =K log 2 +2,-2. (4) 

x 
Since this is a transcendental equation which cannot be solved in 


a general closed form, we shall consider two special cases. 


Case 1. K >>x. In this case equations (2) and (4) give the fol- 
lowing result: 


2 = Ar, —————_, a#)}; a=Ag,te-™*, ag =b. (5) 
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We shall suppose that when z exceeds a threshold value 2* an ef- 
fect, such as a lethal effect, occurs. When it does occur, it will 
take place at a time ¢, which depends on the dose @o- The thresh- 
old dose, x*, is that value of @, for which the maximal value of z 
equals 2*. The quantity x* is generally the most useful single 
number associated with this particular drug and its effect. It is an 
LD50 if the effect is a lethal one. Actually, since we deal with a 
population of individuals, x* represents only a median value, and 
an additional parameter is required to determine the variability of 
the observed effect. The variability is also an important quantity 
but we shall not consider it here. The measured value of x*, the 
LD50, depends to some extent on the time limit which is chosen. 
Practically, however, unless there are appreciable numbers of de- 
layed deaths, the measured value is independent of the time chosen, 
provided the time is long enough. If we find substantial numbers 
of delayed deaths so that the temporal distribution of mortality is 
bimodal, we would suspect that we are dealing with two effects 
which might be separated by observing various pathological changes. 

It follows from (5) that when the dose is increased the median 
time of death should decrease to zero. However, there will always 
be delays which are, to a first approximation, independent of the 
dose. Hence we may add a time ¢, to the value given by (5) to 
obtain the expected death time. The principal temporal aspects of 
the mortality time as a function of dose can be obtained by ob- 
serving when deaths occur at an LD50 dose, at twice that dose 
and at ten times that dose. The parameters can be calculated from 
such data in the following way. Let ¢, be the time when 2 in (5) 
reaches its maximum: 


i 
a0) 


a 
AR a t, =Va,a=b. (6) 


We now calculate z*, the median threshold dose, in terms of the 
parameters for this case. When ¢=¢, from (6) and 2 = 2* in ex- 
pression (5), then by definition z, must be x*. Solving for z* we 
find 


(a — b) 2* ; aez* 
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The value ¢, given by expression (6) is approximately equal to 
the mean time at which z reaches its maximum in the population 
unless a >> b, in which case ¢, becomes large compared with 1/a 
while the mean time of effect will never exceed 1/a by a very large 
amount. 

If ¢,, tg, and ¢,, are the mean times of death corresponding to 
one, two and ten times the LD50 dose (Figure 1) then the following 
relations should hold among the corresponding experimental 


variables: 


Since the fourth column is approximately 1/6, we may write 
t.. = 1.2¢,, —0.2¢, (8) 


from which to estimate ¢,,. Having estimated ¢,,, we may estimate 
a/b from columns two or three. Column three may not give a good 
estimate if ¢,, and ¢, are too similar in value. Having an esti- 
mate of a/b we may obtain a from the relation a = y/(t, —¢,,), the 
values of y being given in the last column. 


DOSE X_/x*—> 


FIGURE 1 
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Case 2. K <<x. We consider next what might be expected if 
the equilibrium constant K is small compared with zx. In this case 
we have 


t= 2) - mE yt (9) 
so that 
Az mE, A 
a= —°(1—e-%) — —2— (at -~14 e7**), (10) 
a a 
and 
at, = log ( + sai (11) 
mE 9 
AmE, [ ax ogee: 
eee ioe {1 |. (12) 
a mE, mE , 


Since these relations are transcendental, we shall only consider 
some special cases in which certain terms can be neglected. In 
applying the results it is necessary to keep in mind the imposed 
inequalities. 

We next give the results for two subcases, 

Case 2a. If az* >> mE, so that ax*/mE, >> log ax*/mE, , we 
find, after elementary calculations, 


- -69 a eee A .105 
bat tea, La Seep aie ei, (13) 
i, =t,, ax* by te aax* 
log log 
mE mE, 
Hence 

i -t 
ee SLO (14) 
me ee 


This value is comparable to the values of column 4 in the above 
table. From (8), or a slightly modified expression, we obtain ¢,, so 
that from (13) ax*/mE, can be calculated. By hypothesis, this 
must be considerably greater than one. Using this value and noting 
that ¢, =¢,-¢,,, the value of a can be calculated from (11). We 
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find a = .69/(t, —¢,,). Similarly from the value of x* the value of 
a/mE , can be calculated and hence also mE,: 


69r* ~.69(t,-t. )/(to-t_) 
EWE ee (15) 


mE = 
imag 
Case 2b. K << x* << mE,/a. If « is small compared with mE ,/a 
but still much larger than Kk, then 


ax ax 
at, = log (1 - se} aie ) (16) 
and 
2 
Be) fosed 9 Ps (17) 
* 28 \nE; 


In this case z is no longer proportional to z,. If twice or ten times 
the threshold dose is given then we find, after elementary 
calculations, 


t, -—¢ £ -—t 
en | oe 10 2S 3 0.05, (18) 
ty genoa = co 
so that 
t,,.-—¢6 
gi ioe Ree (19) 
t,-t,, 


This value is only slightly larger than the values in column four of 
the ahove table. In this case the ratios in (18) must have the 
numerical values required. If so, then from (16) we may estimate 
t)/mE, since t, =t, —t, when 2, =a*. But since z* is known, 
mk, is estimated as x*/t_. The value of a cannot be found though 
it must be much less than mE, /x*. Note that the values in (18) are 
not much different from Case 1 in which 


7 = 1/t or mE, = K/t« (20) 


Il. Effect of two equal doses given T units of time apart 
We consider next the effect of giving two equal doses of the drug 
X, T units of time apart. Let each dose be 4 Fx* so that F is the 


dose in units of 5 a*, We may find, for each value of 7, the value 
of F which is necessary for a threshold effect. It is evident that 
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1<F <2. If, for example, F = 1.5 for a certain 7 then we may say 
that the total spaced dose is only = as effective as a single dose, 
anes phe 

because the total dose (1. 5 + 1.5 eate a Bay a is one and a half 
times as large as the corresponding single dose. Thus F7! is the 
relative effectiveness which we wish to calculate as a function of 
the time interval 7. 

Case 1. K>>x. In this case we first calculate z and z for the 
time ¢ = T whenz, = 4 F2*: 


x(T) = 5 Fx* She ee (21) 
a(T)= 4 Ha tee ee (22) 
a-b 


Since for ¢> T, x = 2(T)e~°', where ¢ is now measured from the 
injection of the second dose, we find 


Ax*F 
zg : ; fers = ene e412 pre CPs eno ye a#b; 
B=2AFe* [Te + (lse tle; a=d. 
The maximum value of z occurs att =¢, given hy 
1 a{l 2") ah 
= og ———__—_ a ; 
Sy ay? Rice fi 
(24) 
: bce wo ale 2% ; 
en C= e 
r a(1+e 97) 


The maximum value z_ of z is now obtained by substituting the 
value ¢ for ¢ in the expression for z. Performing these substitu- 
tions and introducing the value of x* from (7) we find 

b a 


Fu! =i * e eT) a-b (4 i e7 oT) a-b a#zb; 


(25) 


aPe=s 
x (1+e"°") exp (Bese t le a=b. 
In Figure (2) values of F~’ are graphed as a function of 67 for 
various values of the ratio a/b. Since we are interested in the 
shape of the curve for the estimation of the ratio a/b, the curves 
in the figure have been made to coincide at F-! = 3/4 by multiply- 
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RELATIVE EFFECTIVENESS F7 


TIME INTERVAL BETWEEN DOSES bT —> 
FIGURE 2 


ing the abscissas by the following factors before plotting: For a/d 
having the respective values 1, 2, 5, 10, and ~ the corresponding 
factors were 1, 1.40, 1.82, 2:0, and 2.5. From the figure it can be 
seen that the ratio a/b can be estimated only very roughly from 
such data in this case. But the values of a/d and especially a 
should be compatible with that obtained for Case 1 of Section I. 

Case 2. K<<x. In this case, after the second injection, z is 
given by the following expression, provided that T < x*/mE,: 


a(t)=22,-mE,T -mE,t. (26) 


From equations (2) and (26), 2(¢) is easily obtained. The time 
t,, at which 2 is a maximum can then be calculated. Introducing 
this expression for ¢,, into the equation 2(¢_) = 2* we obtain a re- 
lation between F and 2*. By hypothesis #, = 4Fa* and 2* = 
2(¢_) when 7 =0. In this way we find the following relation be- 
tween F and z*: 


w’(F —1) = log (2 + Fa’ + Fa’e*) - log 2(1 + 2x’), 

27) © 
a’ = az*/mE, - (2) 
It can be shown that for x” >> 1 and for T < a*/mE,, the above re- 
duces to 


a E 
Pet Pope (28) 
x a* 
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where f=1. But if F~' is greater than about 0.6, the expression 
(28) is quite satisfactory for any value of x’ if we introduce cer- 
tain numerical values for f. When zx’ has the values 050, 0515-0> 
100, the corresponding values of f are 0.5, 0.61, 0.63, 0.84, 1.0. 
If the time scale for the case in which f =1 is multiplied by a 
factor of 5.2, then the function F~' is very nearly the same nu- 
merically as that of Case 1. For comparison with Case 1, the 
quantity F~' for this special case is included in Figure 2 (dotted 
curve). Note that in this case the abscissa should be (a/z’)T in- 
stead of 67. It can be seen that from data in which T is varied we 
could not easily differentiate Case 2 from Case 1 with a>> b. If 
values of F~' are desired for T > «*/mE, a new expression must 
be derived. But even for xj as low as 0.1, F~' has dropped to 0.61 
before we need a new relation. It may be noted that, when 2’ is 
about one, F~* is actually lower than shown by the dotted curve. 
When 2’ = ~, F~! is about 0.57 instead of 0.62. The curves for 
2’ equal to 0.1 and 10 diverge considerably less. 

Except for the uncertainty in f, we can obtain mE ,/z* from (28) 
with known values F and JT. Hence multiplying by the median ef- 
fective dose x* we obtain an estimate of mE,. By comparing this 
value with estimates of mE, from Case 2a or 2b of Section I, we 
may see whether the data are compatible with this case. 


III. Constant dose rate, R, for K >> x 


The case of a constant dose rate R for Case 1 in which K >@ 
can be treated easily. Two situations arise. In experimental 
situation A, we assume the dose rate is maintained until the effect 
is produced. In situation B, we assume that the dose rate is main- 
tained for a time 7. In either case, we shall solve only for the time 
which gives a median effect. 

Experimental situation A. Dose rate maintained until effect is 
produced. As long as the dose rate is maintained we have 

ea PS 3 wher ATs (29) 

dt K 
We may obtain an expression for x from equation (29). Substituting 
this expression for z in (2) we obtain a differential equation which 
can be solved for 2. If we set z = 2* and ¢ = ¢, we have a relation 
between the dose rate R and the time ¢, which is required for 2 to 
reach its threshold value. The time ¢, is the median survival time 
less ¢,,. If the dose rate is too small, 2 will never reach 2*. If 


-: 
ae 
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we let R* be the least value of R which permits z to reach z*, we 
find that the relation between FR and ¢, can then be written in the 
following form: 


R* be “*s — ae °'s 
Frou 14 == ypu oa (30) 
a _ 


From this expression we can graph R/F* in terms of ¢, for various 
values of a and 6, then invert, giving ¢, as a function of R/R*. 
If we are interested in the accumulated dose as a function of dose 
rate or of total time, this can be obtained as shown in Figure 3. 
Note, however, that these values do not include R¢,, and the time 
t, does not include ¢,,. In spite of this the accumulated dose at 
first decreases but then increases with the time ¢, . 

Experimental Situation B. Dose rate maintained for a time r. In 
this case we obtain 2(r) from (29) and then 2(r) from (2), with 
2(0)=0. After the time 7, x(t) is given by (29) and 2(Z) by (2), 
with the initial condition z = x(r). Solving for the time ¢, at which 
z 1S a maximum, introducing this value into the expression for 


TIME bT 
---- TIME bf, 


nN 


ACCUMULATED DOSE 


TIME 
FIGURE 3 
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z(¢>7),:and setting this expression equal to z2*, we find the fol- 
lowing relation between 7 and R: 


R* a cue 
ee Ba (1 cf eo ee (1 ae et ES lea ; (31) 
R 
Here again we can obtain 7 in terms R/R* for various values of 
a and 6. Cumulative dose curves are shown as a function of total 
time 7 in Figure 8. The accumulated dose always increases with 
time r. 


IV. Median threshold dose for two drugs having a common effect 
and a common detoxification pathway 


If two drugs X and Y, which are injected either simultaneously 
or a known time interval apart, are inactivated through a common 
enzyme then we may consider the following reactions: 


Sk EX) Bak, 
Wo ste Y) > bay, 


Case A. Drug Y detoxified slowly compared with X, m’ << m, 
K’ << K. We shall assume for simplicity that m’ <<m. If K >> Kk’, 
the concentration y of Y is much smaller than that of (FY), and we 
have 


a being defined by (32). Let &z* be the amount of X injected and 
ny*—the amount of Y so that € and 7 are the amounts relative to 
the respective median lethal doses, but such that the combined 
dose results in a median effect. Because we have assumed that m 
is small, therefore, y* = az*/A’, where A’ measures the effective- 
ness of Y. The effect of drug X is most easily considered sepa- 
rately for the cases K >> « and K <<z. When z has values com- 
parable with K the equations are much more complex. 

Case Al. K >>x. Integrating (32) and introducing x(t) into the 
following expression, which is similar to (2) but in which the ef- 
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fect depends upon both X and Y, we find 2(¢): 


2(t) = eee e* (Ar +A’y)dt. (33) 
0 


When &a* and ny* are injected simultaneously the maximum value 
of z is by definition equal to 2*; therefore, we find the following 
equations for the determination of the relation between € and 7: 

t -a 


a-—b ae? ‘my & 


a-a brs 


at 


oe -(l1-e ™)n. (34) 


aca 35 
| 65) 


(36) 


and 


b y* az* 
a = ~ ha 
1+ Bn | ee Ae 
In Figure 4 are shown some curves relating € and 7. Also shown 
are some curves in which there is a time T between injections, Y 
being given first. The effect of the interval T will be discussed 
again later. The straight line from (0,1) to (1,0) corresponds to the 
injection of two drugs having identical effects and having the same 
time constants, such as would occur with two simultaneous in- 
jections of the same drug. As a becomes considerably less than } 
for the case in which the effect of Y on X is negligible so that 
B = 0, the curve approaches a rectangle. This corresponds to the 
situation in which the processes are independent (cf. independent 
joint action, Bliss, 1939), even though by hypothesis the two drugs 
act through the same final path. The rectangular curve will, how- 
ever, only be approached if correction is made for mutually exclu- 
Sive probabilities so that at the point (1,1) the expected effect will 
occur with a probability of 3/4. 
In this situation we may say that Y potentiates the effect of X 
because of the conditions imposed on the parameters. Note, how- 
ever, that the ratio of the threshold doses, x*/y*, is equal to the 
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b=Sa, A- 0 T=0 
S 


of 3 4 6 8 1.0 
§é—> 


FIGURE 4 


ratio A’/A multiplied by a number lying between 1 and e. The 
number is one when a >> } and e when a= 5b. All the curved lines 
of Figure 4 approach the €~axis at smaller angles than they do the 
n-axis, if they reach the latter at all. 

We may refer to the effect of two drugs as being directly addi- 
tive under given conditions when the result is represented by a 
point (€,7) lying on the straight line +7 =1. If the point repre- 
senting the interaction lies to the left or below this line we may 
refer to the interaction as being more than additive. If it is to the 
right or above, we have less than additive potentiation. If the 
point lies to the right or above the rectangle (0,0), (0,1), (1,1), 
(1,0), this may be called negative potentiation. It will be discussed 
subsequently. It can be seen from the figure that in our model the 
resulting effect may be greater or less than additive, depending on 
the time 7 between injections. Similarly, for a given time, the 
effect may be more or less than additive depending on the ratio of 
the drugs used. 

Case A2. K << x. In this case we find that € and 7 are related 
by the following expressions: For azj >> mE) , 


at 


(ce ™-e ™E+(1-e ™)n, (37) 


a-@ 
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where 


log ~ ( + ee (38) 


a* is given by (7), and y* = az*/A’. 
For ax* << mE, we find the following transcendental expression 
relating € and 7: 


Ta — 'O 


2AmE -o —at —at 
1s je le ere ee eee (39) 


a? 2* 
t,, being given by (38) with at? = 22*mE,/A. The above expres- 
sions do not hold when 7 is too small. 

Case B. Special case in which both drugs have comparable rates 
of detoxification; K, K’ both small; m’K = mK’. In this case we 
find from (82) and the corresponding expression for y: 

og Cae mm°E yt 


2S Se 
Yo Lo MYy +M Lo 


(40) 


The maximum value of z is given by setting the derivative of 2 in 
(33) equal to zero and solving for the maximum z_: 


20 = ato + Ao - ye AtntiA ad + A’yo) log (1+ ae pee (41) 
- te a To Yo mE, mE, 
a | —— + 
mE, mE, 


Case B1. If ax,/mE, >> 1, ayy/m’E, >> 1, then we find €+ 7 =1 
so that the potentiation corresponds to uniform additivity, that is, 
the drugs add in all proportions. 

Case B2. If axy/mE, << 1, ay,/m’E, >> 1 we find € and 7 re- 
lated by 


2mE,A 
This is a straight line starting at (0,1) with a steep slope and 
intersecting the abscissa far to the left of one. But, as pointed 
out above, when 7 becomes small, the approximation does not hold. 
Hence, the relation between € and 7 in this case resembles that 
shown in Figure 4 for Case Al in which 8B is large. 
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Case B3. If ary << mE, and ay, << m’E, we find the following 


relation 
( mA, la mA’ } 
Tap MN rare on) Raia ey 


In this case the curves are symmetric about the line +7. Only 
One parameter is involved. Expression (48) is illustrated graphi- 
cally in Figure 5 for several values of the parameter, mA/m’A’. 
When the parameter has the value one, the drugs add in all 
proportions. 


V. Potentiation in the case where detoxification may be 
neglected and where both drugs may be bound 
to a common substance 


Let K and Kk’ denote the eouilibrium constants for the combina- 
tion of X and Y respectively with some substance. This substance 
may be an enzyme or it may be a substance which combines with 
X or Y in the presence of some enzyme which is plentiful and which 
acts rapidly. Let both X and Y have great affinity for the sub- 
stance. In our notation this means that K and Kk’ are small. Let 
E, be the total amount of the substance. Assume that the equilib- 


Cae 


FIGURE 5 
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rium concentrations are rapidly reached and neglect detoxification 
other than that which results from X and Y being bound to the sub- 
stance. In this case we find the following equation giving the rela- 
tion between &é and7, if K =k’, 


(14 2°40 E+ (14 2’)n¢ 2) [(1t 2° +e E+ (i+ (1+ 2’)n-(1+ Ol = 
(1+6)2’, (44) 


where «=(4—A’)/A, 2’ =az2*/A’E,. It can be shown that if 
«= 0, +7 =1 and the result is uniform additivity. Furthermore if 
(1 + «)2’ = 1, the result tends to be nearly the same as if e=0. 
If (1 + ©) and 2’ are much greater or much less than one, the result 
is ‘‘greater than additive.”” In particular if (1 + «) = 2’, then we 
find that the relation between 7 and € is the same function of 
2’=1+e as it is of 1/2’. The value of &€ when &=7, 1.e., if 
equally toxic doses are used, has the following value: 


1 Gea 
413 4 ADate ae = 


E=7n for 2° =1+e, (45) 
so that € lies between 1/3 and 1/2. However, as 2” becomes large, 
the relation between 7 and & becomes very asymmetric so that the 
point (+, 0) is nearly on the curve. 

If we do not restrict ourselves to 2” = 1 + « then we have 


24+ 8¢4¢7 4 292 + 27 + 42 + Seg 


Te Ee (2+22+6)(2+22+2+e2) — Ce) 


There are now two conditions under which € = 7 — 0 so that 
the potentiation becomes very large. They are «>> 2’? >> 1 and 
(1 + &) << 2’? << 1. These conditions, which favor great potentia- 
tion, can also be written as (A/A’) << (az*/A’E,)* <<1 and 
(A/A’) >> (az*/A’E,)* >> 1, so that the intrinsic toxicity of the two 
drugs differs greatly. 

It may be noted that the points €=0, 7 = 1 and €=1, 7 =0 are 
always points on the 7(&) curve when K = K’, although the curve 
may under some conditions approach the &-axis at a point con- 
siderably less than one. 

If K =k’ and if (K — K’)/K =6 then € and 7 are related by the 
following expression: 


[e+ d)(L+e+ 2’) &+ (e+ d5)\(1 + 2’)n - (eo? + 82° 4 5 4 €d)] x 
[e+ dL +e + 2&4 (1+ (e+ d)(14+ 2’)n -(14+ (e454 S2’)] = 
(1+ ee? 2°(1-8). (47) 
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If both K and K’ are large then since neither X nor Y is bound to 
an appreciable extent, it is to be expected that €+=1. If, for 
example, K is very large, then we find »(€) as follows: 


(K° + 2° +E, -2°€) 


Kit eal & 


n=(1-€) (48) 


Since the right-hand side is greater than (1 - €), the result is that 
the potentiation is uniformly less than additive. 


VI. Conditions for additivity in all proportions 


In view of the results given above, we would not expect uniform 
additivity, i.e., additivity in any proportion of drugs administered, 
except under special conditions, even though by hypothesis we are 
considering a single common effect. However, there are a number 
of conditions which give rise to uniform additivity. Some of these 
are given in the following list: 

A. Common effect and common detoxification paths. 

1) K = K’, A= A’*;m =m’. 
2) Km’ = K’m, K and Kk’ small, x and y5, >> 1. 
3) K/K’ =m/m’ = A/A’, K and kK’ small, x6 and yf << 1. 

B. Common effect, negligible detoxification, but temporary in- 

activation by chemical combination. 
1) K and Kk’ large. 
2) K and K’ small, A =A’. 


VII. Effect of time interval on potentiation 


In the above cases in which the time constants were not all 
neglected, it is to be expected that the order of injection and the 
time interval involved will influence the outcome. In a few cases 
calculations have been carried out to illustrate some of the ex- 
pected relationships. These are shown in Figure 6 for the Case Al 
of Section IV in which equal relative doses of the two drugs are 
administered. It can be seen that the order is important and that 
potentiation may be less than additive for simultaneous injections, 
while greater than additive if the proper drug is given a suitable 
time before the other. If the detoxification of Y is not completely 
neglected, the curves reach a maximum and decrease toward one 
half for long intervals (dotted part of curve, Fig. 6). 


SA 
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VIII. Antagonistic interaction 


Next we consider a case of two drugs, each producing a lethal 
effect through a different process. Let X act through Z, just as 
discussed. Similarly let Y act through W,:w being the concentra- 
tion of W. We shall illustrate the theoretical results for three kinds 
of interaction, but in all cases we shall neglect the inactivation 
of X and Y in order to simplify matters. In the first case we shall 
Suppose that the drugs interact directly. In the second we con- 
sider the situation in which the quantities Z and W interact. In 
the third we suppose that X and W interact at the same time as do 
Y and Z. 

Case 1. Direct interaction. Let the effect of the interaction be 
such that X and Y are eliminated by the formation of (XY), the 
equilibrium constant of the reaction being K,- Let X act through Z 
where 2 is given by (2) while Y acts through W, w being given by 
an expression similar to (2) in which @ and A are replaced by a@ 
and A. Then a =, ~(XY) and y = Y y ~(XY). Furthermore, in this 
case, if w* is the threshold lethal value of w, 


a* = az*/A, y* = aw*/A. (49) 
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Suppose that a certain quantity 2, of X and Y> of Y, when given 
together, produce a median lethal response z* as a result of X. 
Let 2, = €’x* and ¥,=n’y*. This defines €’ in terms of 7%. 
(Loewe, Joc. cit.; Gaddum, 1948). Using (2) together with the re- 
lations just discussed and setting 2 = 2* we find a relation be- 
tween €’ and 7’, which is represented by a straight line passing 
through the point €* = 1, n’ = 0. This illustrated in Figure 7. 

In a similar way we can find a relation between &”’ and 7’’, con- 
necting the individual doses when the drugs are given in combina- 
tion but where the lethal threshold is due to W. This is also repre- 
sented by a straight line, but it passes through the point €’’ = 0, 
n’’ = 1. The two lines intersect at a point (€,,7,), the coordinates 
€, and 7, of which are given by 


Pen ie y*/Ke n, =1+2*/K,. (50) 


It is clear that only the part to the left of and below these two 
lines is of interest. Thus the desired relation 7(&) is made up of 
the parts of n’(&’) and 7n’’*(€’*) between (€,,7,) and the axes. 
Note that at the point (€,,7,) the mortality, uncorrected, would 
be 75% if the sensitivities of the two effects are uncorrelated. In 
this case it is seen that we can calculate K,, the only parameter, 
from experimental values of €, and y* by using the first expression 
of (50). However, K, can also be calculated from the second ex- 
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pression of (50). If the results are compatible with this case, the 
data fall along approximately straight lines, except near (E57: 
unless corrected. Note that for this case 


ean (51) 


the ratio of the median lethal doses. Hence the distance that the 
point (€,,,) is from the origin depends upon the parameter K, - 

Case 2. Interaction of the intermediates Z and W. In this case 
we consider a situation in which the intermediates are eliminated, 
at least temporarily, by a reaction Z +W >» (ZW). In this case it 
is simpler if we use the following differential equations for 2 and 
w: 


d 

oe Agr -az-—nzew, (52) 
dt 

dis v= 

— = Ay-GQw-nzw. 53 
7 y (53) 


Since we are only interested in the quasi-steady state we fol- 
low a procedure similar to that used in Case 1. The relation be- 
tween €’ and n’ can be obtained as follows. Let the time deriva- 
tives in €’ and 7’ equal zero, e = 2, = €’x*, y=y, =n’ y* and set 
z2=z2*. Eliminating w and introducing (49) gives an expression in 
which 7’ is linear in €’ and which passes through the point €’ = 1, 
n’ = 0. Similarly, letting €’’, n’’ represent the case in which the 
threshold w* is reached, we let x = €’’x* and y=7n’’y*. Setting 
w = w* and eliminating z we find an expression in which 7’’ is 
linear in €°’ and which passes through the point €” = 0, n’ = 1. It 
can be shown that the lines n’(€’) and n’’*(€é’’) intersect as shown 
in Figure 7 at a point (€,, 7.) given by 


nAy* nAa* 
€,=1+ ae" Ng =1lt+ . 


aa (54) 


aa 
In this case the expected results differ from Case 2 only in that 
two parameters are now involved. For parameters we conveniently 
take nA/aa@ and A/A. The latter is the ratio of (n, — 1) to (€, - 1), 
times the ratio of the median lethal doses. It does not have to 
equal one in this case as it does in Case 1. 
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Case 3. Each drug acts on an intermediate. If we suppose that 
the interaction can be represented by 


Mate Wre a (XW) eo V4 % 2s: (YZ)', 


we find from the method used in Case 1 that 7(€) is again made up 
of two straight line segments starting respectively at the points 
(0,1) and (1,0) and intersecting at a point (€3,73) given by 

Ay* An* 


=e | Ng =1i+ : 
ak, ak, 


at (55) 
Here again two parameters are involved and the result is identical 
with that of Case 2. 

It is possible that a lethal effect can be produced by more than 
one process. If the threshold dose for process one is substantially 
lower than the threshold for process two, which results from the 
Same drug, the latter can be ignored for practical purposes when 
considering the drug alone. However, if a second drug in any way 
reduces mortality that would result from the first process of the 
first drug, then the second process, if unaffected by Y, may now 
produce its lethal effect. The result of such a situation is also 
illustrated graphically in Figure 7. Here it is assumed that Y also 
has a second lethal effect requiring a median dose F’y* and that 
this effect is not influenced by X. In this case the ‘‘corrected’’ 
n(&€) curve would be a polygon (0,1), P, Q, (1,0), provided that any 
second lethal effect due to X would have a median threshold 
greater than ¢,2*. 


Discussion 


Define the potentiation ratio P as the distance from the origin to 
a given point ¢,7 divided into the distance along a line from the 
origin through the point €,7 to the line + 7 = 1. The ratio is then 
the dose required on the assumption of uniform additivity divided 
by the amount actually required to produce a threshold effect. If 
this ratio is plotted against the proportion of either drug, e.g. 
€/(€-—7n), both drugs being measured relative to their effective- 
ness, it is found that the resulting curves, which pass through the 
points (0,1) and (1,1), will be straight lines for uniform additivity. 
The curves will be concave downward for simple cases of greater 
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than additive potentiation, and concave upward for cases in which 
there is less than additive potentiation or antagonistic action. It 
may have an inflection point so that the ratio is greater or less 
than one depending upon the proportion chosen as well as upon the 
interval between administrations of the drugs. Thus one cannot 
speak of the value of the potentiation ratio without specifying the 
relative doses and the time interval used. 

Under some conditions the largest value of P may be a useful 
quantity to know for a given pair of drugs. If the potentiation 
ratios for the various cases discussed above are plotted as func- 
tions of €/(€ +7) then it is found that the maximum value will 
rarely be much more than 10% of the largest value at € to 7 ratios 
of 1:4, 1:1 or 4:1. The exception to this occurs where the series 
of P values steadily increase or decrease, the largest value repre- 
senting a potentiation ratio greater than about five. In this case 
the largest value of P may be as much as 20% less than the the- 
oretical maximum. Hence we would anticipate that the use of a 
pair of drugs in the proportions 1:4, 1:1, and 4:1, the doses being 
measured in units of relative effectiveness, would enable a fairly 
satisfactory estimation of the maximum potentiation ratio for a given 
set of conditions including the interval of time between doses. 

The problem of variability has only been referred to above. If 
two drugs have actions identical in all respects, then the regres- 
sion coefficients of the log-dose-response curves will also be 
equal. However, if the two drugs have the same final effect but 
with one showing a predominantly linear detoxification resulting 
from saturation of the detoxification enzyme, while the second 
shows an exponential detoxification in the toxic range (cf. Section 
IV, Case B2), then the regression coefficient of one will be twice 
that of the other. When both drugs are used the regression coef- 
ficient may take on intermediate values, even though both drugs 
act through the same mechanism. 


The author is indebted to Dr. K. P. DuBois for valuable 
suggestions. 

This work was aided by a grant from the Dr. Wallace C. and 
Clara A. Abbott Memorial Fund of The University of Chicago. 
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MINERALISCHER NAHRSTOFFE 
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Fur die Praxis der Pflanzenernahrung ist es wichtig, zu wissen, in 
welcher Weise die Ertragsbildung von der Konzentration eines minerali- 
schen Nahrstoffes in der Umgebung der Pflanze abhangt. Da nur diejenigen 
Nahrstoffmengen das physiologische Geschehen in der Pflanze unmittelbar 
zu beeinflussen vermogen, die sich in der Pflanze befinden, wird angenom- 
men, dass das Wachstum zum Zeitpunkt ¢, d.h. die Geschwindigkeit der 
Trockensubstanzzunahme zu diesem Zeitpunkt, eine Funktion der zur 
Zeit tinder Pflanze enthaltenen Nahrstoffmenge ist. Diese Nahrstoffmenge 
wird naturlich im Intervall vor dem Zeitpunkt ¢ aufgenommen. Deshalb 
und auch noch aus anderen Gruinden hangt das Wachstum zur Zeit ¢ davon 
ab, wie die in der Umgebung der Pflanze herrschende Konzentration des 
betrachteten Nahrstoffes in demjenigen Zeitintervall verlauft, das sich 
von der Aussaat bis zum Zeitpunkt ¢ erstreckt. Die angegebene Annahme 
fuhrt zusammen mit einigen weiteren naheliegenden Annahmen zu einem 
Ansatz, der Ergebnisse liefert, die in verschiedener Hinsicht gut mit der 
Erfahrung tibereinstimmen. Jedoch gibt es auch noch Widerspruche 
zwischen Theorie und Erfahrung. Durch weitere Ausgestaltung der Theorie 
lassen sich diese Widerspriiche beseitigen. Es wird angeregt, Versuche 
durchzufiihren, deren Resultate Hinweise fiir die weitere Ausgestaltung 
der Theorie liefern. 


In einer friiheren Mitteilung (Kaltofen, 1957) wurde von der na- 
heliegenden Annahme ausgegangen, dass das Wachstum dy/d¢, d.h. 
der Differentialquotient der Pflanzentrockenmasse y nach der Zeit 
t, der Trockenmasse zum Zeitpunkt ¢ proportional ist und ferner 
davon abhangt, wie hoch der relative Gehalt der Pflanze an dem 
gerade betrachteten mineralischen Nahrstoff zum Zeitpunkt ¢ ist. 
Unter relativem Gehalt an dem betrachteten Nahrstoff ist der durch 
die Trockenmasse dividierte absolute Gehalt der Pflanze an dem 
betrachteten N&hrstoff zu verstehen; der relative Gehalt werde mit 
»” bezeichnet. 
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Die angefihrte Annahme fiihrt zu dem Ansatz 


dy 

ap YH (2): (1) 
Dieser Ansatz gilt zundchst nur fiir das Intervall ¢, <¢ <¢,, in dem 
einerseits das Saatkorn praktisch keine Rolle mehr spielt und 
andererseits Reife- und Absterbevorginge noch von untergeord- 
neter Bedeutung sind. Es ist also noch zu untersuchen, welche 
Form ein Ansatz hat, der im ganzen Vegetationsintervall zwischen 
der Aussaatzeit ¢= 0 und dem Zeitpunkt des vélligen Absterbens 
der Pflanze ¢=¢, gilt 

Zunaichst kann man die Giiltigkeit der Gleichung (1) auf das 
Intervall 0 <¢<z¢, ausdehnen, indem man festsetzt, dass y zu 
Beginn der Vegetationszeit nicht die gesamte Trockenmasse des 
Saatkornes darstellt, sondern nur die Trockenmasse des Keimlings. 
Dementsprechend ist unter » zu Beginn der Vegetationszeit der 
relative Gehalt des Keimlings an dem betrachteten Nahrstoff zu 
verstehen. Auf Grund dieser Festsetzungen zahlen alle Teile des 
Saatkornes, die nicht zum Keimling gehéren, wie Endosperm usw., 
zur Umwelt der Pflanze. Pflanzenarten, bei denen die Keimblatter 
zu Speicherorganen ausgebildet sind, seien vorlaufig ausser Betracht 
gelassen. 

Will man des weiteren einen Ansatz gewinnen, der auch fiir 
t>t, gilt, so sind die Erscheinungen des Reifens und Absterbens 
zu berticksichtigen. Im folgenden seien nur einjahrige Pflanzen 
betrachtet. Bei vielen einjahrigen Pflanzen sterben schon Blatter 
ab, lange bevor die Pflan¥e zur Reife gelangt. Die zeitliche Reihen- 
folge, in der die einzelnen Blatter absterben, entspricht dabei 
meist der Reihenfolge, in der die Blatter gebildet wurden. Zuerst 
sterben die zuerst gebildeten Blatter ab, usw. Aber auch wenn 
keine Blatter schon friihzeitig absterben, wird die ‘‘Vitalitat’’ der 
Pflanze mit Fortschreiten der Vegetationszeit immer geringer. Das 
Wachstum je Einheit Pflanzenmasse nimmt im ganzen Intervall 
0 <¢<¢, monoton ab (Rippel, 1925), bis es schliesslich zur Zeit 
der Totreife den Wert Null erreicht. 

Den angefiihrten Gegebenheiten kann man dadurch Rechnung 
tragen, dass man y in (1) mit einem variablen Faktor A = (2) 
multipliziert. Man erhalt dann die Gleichung 


dy 
a ye). (2) 
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X ist eine monoton fallende Funktion von ¢ mit A(¢,) = 0. Weiterhin 
kann man \(0) = 1 setzen. Mit (2) ist ein Ansatz gewonnen, der im 
ganzen Vegetationsintervall 0 <¢< t, gilt. 

Fiir die Praxis der Pflanzenernéhrung ist es nun wichtig, zu 
wissen, in welcher Weise v von der Konzentration des betrachteten 
mineralischen Nahrstoffes in der Umgebung der Pflanze abhingt. 
Denn wenn man dies weiss, kann man wegen (2) auch Aussagen 
dariiber erstatten, in welcher Weise das Wachstum dy/dt bzw. die 
Trockenmasse y von der Konzentration des betrachteten Nahrstoffes 
in der Umgebung der Pflanze abh&ngt. 

Ankniipfend an friihere Uberlegungen (Kaltofen, 1957), erhalt 
man fiir » den Ausdruck 


oe == (NO) “W() oe) at | 3 
rar +| y) du) dt | . (3) 


Hierbei ist N = N (2) der absolute Gehalt der Pflanze an dem betrach- 
teten mineralischen Nahrstoff. W(y)¢(u) ist die Geschwindigkeit, 
mit der die Pflanze diesen N&hrstoff aus der Umgebung aufnimmt. 
W ist die Oberflache, die den betrachteten Nahrstoff resorbiert; W 
ist von y nicht unabhadngig. Die Argumentfunktion wu = u(t) be- 
schreibt die in der Umgebung der Pflanze herrschende Konzentra- 
tion des betrachteten Nahrstoffes, die beliebig verlaufen kann. 
¢(u) ist eine passend gewahlte Funktion von w. 

Im folgenden sei an Hand einiger Beispiele gezeigt, dass der 
Ansatz (2) in Verbindung mit (3) im Prinzip brauchbar ist. 

Da iiber die Funktionen d(¢), / (v), W(y), (wu) noch nichts Ge- 
naues bekannt ist und um die Berechnungen vorerst méglichst 
einfach zu gestalten, sei vorlaufig A(¢)=1-¢t,~'t, (v) = av, 
W (y) = by, d(u) = cu gesetzt, wobei a, b, c positive Konstanten 
sind. Damit erhalt man aus (2) unter Beachtung von (3) die Integro- 
Differentialgleichung 


ad t 
= a(1~ ty”) [No + bo | uydt|. (4) 


Selbstverstindlich stimmen die Funktionen y = av und ¢ = cw nur 
dann einigermassen mit der Wirklichkeit tiberein, wenn der Wertbereich 


der jeweiligen unabhangigen Veranderlichen, v bzw. u, hinreichend 
eingeschrankt ist. Die Funktion ~(v) beispielsweise steigt zwar 
mit wachsendem »v zundchst an; von einem bestimmten Wert v an 
fallt sie jedoch wieder ab. 
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Mittels Iterationsverfahren kann man Naherungslésungen der 
Gleichung (4) berechnen. Zundchst sei nur der Fall betrachtet, 
dass u(t) im ganzen Intervall 0<¢<¢, identisch konstant ist, 
wobei der triviale Spezialfall, dass w(¢) im ganzen Intervall0 < ¢<¢, 
gleich Null ist, ausser Betracht gelassen werden kann. 

Setzt man fiir y in (4) den gegebenen Anfangswert y(0) ein und 
integriert man daraufhin die Gleichung (4), so erhaélt man unter 
Einsetzung des Anfangswertes y(0) als Integrationskonstante die 
Néherungslésung 

y=A4,+A,t+A,t? +A,e? (5) 
mit den Koeffizienten 


A. = y (0), 
A, =aN(0), 


A, = 4 abouy(0)- 4 aN(0)tg', 


2 2 


A,=- 7 abeuy(0)t,-'. 


Diese Naherungslésung ist jedoch noch zu ungenau. Es empfiehlt 
sich daher, die nachstfolgende Iteration zu berechnen. 

Setzt man in (4) fiir y die N&éherungslésung (5) ein und integriert 
man dann die Gleichung (4), so erhalt man, wenn man wieder y(0) 
als Integrationskonstante einsetzt, die Naherungslésung 


y=B,+B,t+...+B,t® (8) 


mit den Koeffizienten 


Ba = 5 abeuy(0)- 4 aM(0)t,-', 

B,=% @bdcuN(0)- 4 adcuy(0)t,-', 

B= ah a2? o2v? y(0) - $a? BouM Oey 

By = gy a@bcuN(O)tp-? - + a2d3 oF ur y(0)t,', 
By = qq 275207 u2 y(0)t,77. 


Die N&herungslésung (6) ist schon wesentlich genauer als die 
N&herungslésung (5). 


Noch bessere Naherungslésungen als (6) erhalt man, wenn man 
das Iterationsverfahren weiter fortsetzt. Die Rechenarbeit nimmt 


WIRKUNG MINERALISCHER NAHRSTOFFE 29 


dann allerdings betrichtliche Ausmasse an. Aus diesem Grunde sei 
das Verfahren abgebrochen. Eine genaue, quantitative Uberpriifung 
der Theorie ist zur Zeit sowieso noch nicht médglich, weil noch 
keine entsprechenden Versuchsergebnisse vorliegen. Die derzeit 
bereits vorliegenden experimentellen Erfahrungen gestatten lediglich 
eine grobe, qualitative Uberpriifung der Theorie. Fiir eine derartige 
Uberpriifung ist aber die N&herungsldésung (6) genau genug. 
Nunmehr sei der Fall betrachtet, dass w(¢) im Intervall 0 <¢<z, 
identisch konstant und positiv und im Intervall¢, <¢<t, identisch 


gleich Null ist, wobei 0 < ¢, < tp, sei. 
Um fiir diesen Fall eine Naherungslésung zu erhalten, kann man 


im Intervall 0 < ¢ < ¢, das angegebene Iterationsverfahren anwenden. 
Fiir die vorliegenden Zwecke ist die Iteration (6) im Intervall 
0 <¢<¢, hinreichend genau. 
Der weitere Verlauf der Naherungslésung ergibt sich aus der im 

Intervall ¢, <¢<¢, giiltigen Gleichung 

dy & 

Tiara me *t)N (¢,). (7) 
Fiir N erhalt man durch Einsetzen der Iteration (5) in die fiir0 <¢ <7, 
giiltige Beziehung 


N(t) = N(0) + bou [vat 


0 


den Naherungsausdruck 

N(t)=C,+Cyt+...+C,t', (8) 
wobei 

C, =N(0), 
, = Seuy(0), 

4 abcuN(0), 

1 _- 
3 = % 2b? 07 u? (0) - ¢ abcuN(0)t, 
=-75 ab? ou? y(0)ty at 


Seecoen Cy Cy 
tt 


4 


ist. Als Naherungsausdruck y im Intervall ¢, t, ergibt sich 
durch Integration der Gleichung (7) der eee 


y=D,+D,t+D,¢? (9) 
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mit den Koeffizionten 

D, =aN(¢,), 

Di, =-+aNM(t,)t, '; 

D,=B, +B ,t,+:..+B,t,°-D,t,-—D,t,’, 
wobei N(¢,) nach (8) zu berechnen ist. 

Alle in der angegebenen Weise berechneten Na&herungslésungen 
der Gleichung (4) sind im Intervall 0 <¢<z¢, kleiner als die exakte 
Lésung der Gleichung (4), unabhangig davon, wieviel Iterations- 
schritte durchgefiihrt werden. Ferner ist leicht einzusehen, dass 
die nach dem n-ten Iterationsschritt erhaltene Naherungslésung im 
Intervall 0 < ¢ <7¢, kleiner ist als die nach dem (n + 1)-ten Iterations- 
schritt erhaltene Naherungslésung. 

In der Abbildung sind Naherungslésungen der Gleichung (4) 
dargestellt, wobei a= 1, b=1, c=1, y(0)=1, N(0)=1, ¢, = 
gesetzt wurde. Man erhalt die Lésung y,(¢), wenn fir 0<¢<¢ 
u(t)=u,=1 ist, Yo (0s wenn fiir 0<tX<t, u(t)=u, =2 ia 
y,(t), wenn fir O<¢<it, u(t) =2 und far 4 t,<t<t, u(t= 
0 ist. Die Lésungen y, (2) und y,(¢) wurden nach Gleichung (6) 
berechnet, die Lésung y, (¢) wurde nach (6) und (9) berechnet. 

Die erhaltenen Lésungen besitzen einen Wendepunkt, was mit 
der Erfahrung iibereinstimmt. 

Beachtenswert ist, dass y,(¢) im ganzen Intervall 0<¢< tp 
stetig-differenzierbar ist, obwohl die Nahrstoffkonzentration u(t) in 


Yo (t) 


Vat, te 


FIGURE 1 
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diesem Intervall eine Sprungstelle hat. Auf Grund der Gleichung 
(4) bleibt dy/d¢ an Sprungstellen der Funktion u(t) immer stetig, 
was auch in der Wirklichkeit der Fall ist. Ein von B. Baule (1919/ 
20, 1953) angegebener Ansatz ist fiir mineralische Nahrstoffe unter 
anderem deshalb nicht brauchbar, weil nach diesem Ansatz Sprung- 
stellen der Funktion u(t) im allgemeinen auch Sprungstellen des 
Wachstums dy/dt sind (Kaltofen, 1956). 

Von Interesse ist des weiteren, dass y,(t,)>y, (¢,) ist, obwohl 
in beiden Fallen das Integral der Nahrstoffkonzentration wu iiber das 
Intervall 0 <¢<¢, den gleichen Wert hat. Die vorliegenden Erfah- 
rungen (Linser und Pelikan, 1952; Atanasiu, 1953; Selke, 1956) 
lassen durchaus vermuten, dass dieses Ergebnis mit der Wirklichkeit 
tibereinstimmt, sofern das Integral der Konzentration wu iiber das 
Intervall 0 <¢<¢, klein ist. Ist dieses Integral jedoch gross, so 
werden die Verh4ltnisse in der Wirklichkeit gerade umgekehrt liegen. 
Dieser Widerspruch zwischen Theorie und Wirklichkeit ist offenbar 
darauf zurtickzufitihren, dass w%=av die Wirklichkeit auch nicht 
mehr ann&hernd richtig beschreibt, wenn v hohe Werte annimmt. Bei 
Verwendung einer Funktion ys (v), die der Wirklichkeit besser gerecht 
wird als % = av, wird man auf Grund von (2) und (3) auch dann 
qualitativ richtige Aussagen erhalten, wenn das Integral von u 
iber 0 <¢<¢, gross ist. 

Es ware sehr zu begriissen, wenn Versuche durchgefihrt wiirden, 
deren Ergebnisse eine genauere, quantitative Uberpriifung der Theorie 
gestatten oder Hinweise fiir die weitere Ausgestaltung der Theorie 
liefern. Zum Beispiel ware es sehr wichtig, auf Grund geeigneter 
Experimente Genaueres tiber die Funktionen y(v) und ¢(u) zu 
erfahren. 
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The dynamics of cell multiplication and differentiation in tissues in a 
steady state and the kinetics of isotope incorporation into the DNA have 
been theoretically analyzed. Equations have been derived, with the aid 
of which the generation time, the life span, and the distribution or rate of 
death of the cells can be obtained if the tissue is in a steady state, i.e., 
if the number of cells is maintained constant by constant, equal rates of 
cell division and cell death and if the mean DNA content per cell is also 
constant. An equation has also been derived which gives the generation 
time in the case of logarithmic multiplication of cells. 

Two special cases have been analyzed: In Case 1, the isotope is con- 
sidered as being introduced into the metabolic system at zero time only; 
in Case 2, the specific activity of the DNA precursor is considered as 
being maintained constant. The use of the method has been illustrated by 
an example in which the generation time and the mean, the median, and 
the mode life span, as well as the curve of the rate of death of leukocytes 
in a patient with chronic leukemic granulocytic leukemia, have been ob- 
tained from the rate of P3? incorporation into the DNA. 

The merits and the limitations of the method are discussed. 


The rate of incorporation of various radioisotopes into the desoxy- 
ribonucleic acid (DNA) of cells has been frequently used in recent 
years in attempts to estimate the life span of cells (Christensen 
and Ottesen, 1955; Craddock, Perry, and Lawrence, 1956; Hamilton, 
1956, 1957; Kline and Cliffton, 1952a,b; Osgood, Seaman, Tivey, 
and Rigas, 1954; Osgood, Tivey, Davison, Seaman, and Li, 1952; 
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Ottesen, 1954). The methods used to derive these estimates from 
such data vary, and a complete mathematical analysis of the prob- 
lem is lacking. It is the purpose of this paper to present the deri- 
vation of equations that can be used to determine the life span or 
the generation time of cells from the rate of incorporation of an iso- 
tope into the DNA when the number of cells is either in a steady 
state or increases logarithmically. The same equations are appli- 
cable to either tissues or free cells. 


Definitions and assumptions. The term life span as used in this 
paper indicates the interval from the time of formation of a cell by 
cell division to the time of autolysis. The life span of the individ- 
ual cells may be regarded as distributed in some fashion around 
the mean life span of the cell population. 

Generation time is used here with the usual meaning, i.e., to 
signify the time interval between identical phases of two consecu- 
tive mitoses. Individual cell generation times may be regarded as 
distributed around the mean generation time. 

The term steady state pertains only to the number and to the 
mean DNA content of the cells present in the tissue under con- 
sideration and specifically means that they are maintained constant. 

The particular DNA element and its isotope which is used as a 
tracer will be called tracer element and tracer isotope, respectively. 

The mathematical analysis presented in this paper is based on 
the following basic assumptions: 

a) The mean life span and the mean generation time as well as their 
distributions do not change during the experiment. 

b) The number of cells either increases logarithmically, or a steady 
state is maintained by equal, constant rates of cell production and 
cell death. This analysis is thus restricted to specific cases where 
there is reasonable certainty that these assumptions are fulfilled; 
otherwise, it can not and should not be used. 

The following additional assumptions, the validity of which is 
supported by the cited references, have been made: 

c) The mean DNA content per cell is constant (Boivin, Vendrely, 
and Vendrely, 1948; Mirsky and Ris, 1949; Rigas, Duerst, Jump, 
and Osgood, 1956; Vendrely, 1955). 

d) New DNA synthesis occurs only in cells which are preparing to 
divide. The DNA content increases in amount and is subsequently 
divided between the daughter cells (Ris, 1947; Thorell, 1955; 
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Vendrely, 1955). The hypothesis that cell division is accompanied 
by degradation of the preexisting DNA and its replacement by twice 
as much new DNA (Stevens, Daoust, and Leblond, 1953) appears to 
be incorrect in view of the more recent work by the same investigd- 
tors and by others which indicates that the preexisting DNA re- 
mains unaltered while an equal amount of new DNA is synthesized 
prior to cell division (Daoust, Leblond, Nadler, and Enesco, 1956; 
Leblond and Walker, 1956; Nygaard and Rusch, 1955; Takagi, 
Hecht, and Potter, 1956). It is further assumed that the DNA of 
each of the daughter cells is half new and half preexisting DNA 
(Taylor, Woods, and Hughes, 1957). 

e) The isotope incorporation into the DNA is correlated with cell 
division and new DNA synthesis (Brown and Roll, 1955; Hevesy 
and Ottesen, 1943; Howard and Pelc, 1951; Osgood, Li, Tivey, 
Duerst, and Seaman, 1951; Smellie, 1955). 

f) The isotope incorporated into the DNA is not exchangeable with 
the DNA precursors, but it returns to the metabolic pool only after 
cell death (Barton, 1954; Brues, Tracy, and Cohn, 1944; Fresco, 
Bendich, and Russell, 1955; Furst and Brown, 1951; Healy, Simino- 
vitch, Parker, and Graham, 19586; Hershey, 1954; Swick, Koch, and 
Tahara, 1955; Takagi, Hecht, and Potter, 1956). 

g) The DNA of dead cells is not reutilized as intact molecules by 
other cells, but it may return to the metabolic pool following 
degradation. 


Mathematical analysis. The rate of incorporation of an isotope 
into the DNA depends on the turnover rate of the DNA and therefore 
on the turnover rate of the cells (vide supra). Thus the first is a 
function of the latter. We shall derive and express the latter rate 
in terms of the life span L, the generation time /, and the number 
of cells. 


Dynamics of tissues in the steady state. In general, the only 
essential characteristic of a tissue in a steady state is that the 
rate of cell formation be equal to the rate of cell death, which in 
turn may or may not be constant. However, we will consider here 
only the special case in which these two equal rates are constant 
and different from zero. In this case we can distinguish two cate- 
gories of cells: (a) Cells capable of division, which remain con- 
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stant in number and to which a generation time can be assigned but 
not a life span, and (b) cells not capable of division, the number of 
which increases linearly from zero at zero time to a constant num- 
ber at the steady state, when the rate of their death equals the 
rate of their formation. To these cells, which we shall call matur- 
ing differentiating cells, a life span can be assigned but not a 
generation time. Thus if a large number P of cells capable of 
division is present in a tissue which is maintained in a steady 
state, and if each mitosis yields on the average g new cells (normally 
g = 2, but it may be higher because of tripolar or tetrapolar mitoses) 
then the rate of formation » of maturing differentiating cells is 


_ P(g-1) 
~ oo, 


r (1) 
where / is the mean generation time 

The population of maturing differentiating cells in the steady 
state will consist of cells of all ages, i.e., from age 0 (zero) to a 
maximum possible age 7. If the probability of cell survival up to at 
least an age 2 is }(z), then the number of maturing differentiating 
cells which are present at the steady state is 


bi P(g fa. 2 1) T 
M = N da. YS. Th . 
I p (a)dz j |  (x)dx 


Integrating by parts and placing ddé(z#)=-—d[1-¢(z2)], where 
1 — d(z#) is the probability of cell death at an age 2, we obtain 


P(g - 1) F : 
4 = PE=? lo (a ein eat scan}. 


We observe, now, that the first term in the brackets is zero (since 
Lim ¢ (#) = 0) and that the second term is the first moment of the 
x-> 

function 1 — ¢() (Cramer, 1946, p. 174) and, therefore, is the mean 
life span L of the cells. 

Thus we have 


(g -1)L 
P ——. (2) 


The total number T of both types of cells at the steady state is 


(g-1)L+] 
= (3) 


M= 


T=P+M=P 
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The turnover rate of the cells \ can also be expressed as a function 
of M or T by substituting P in equation (1) from equations (2) and 
(3), whereby we obtain 


M 
a L’ (4) 
and 
(g-1)T 
*"@-DLeT 2 


It is apparent that, if any four of these seven variables (T, P, M, 
L, 1, g,; and A) are known, the other three can be determined with 
the aid of these five equations. 

One mechanism by which a steady state can be maintained has 
been proposed by Osgood (1953; 1954; 1957) and is illustrated in 
Figure 1. However, this mechanism is not the only one possible. 
Figure 2 illustrates another mechanism in which the cells capable 
of division fall into two types: (a) the mother cells and (b) the 
multiplying differentiating cells. According to this scheme, when 
a number of mother cells divide they yield an equal number of new 
mother cells and a number of multiplying differentiating cells. The 
latter can subsequently divide only a limited number of times, fol- 
lowing which all of their progeny are converted to maturing differ- 
entiating cells. In this case if m is the number of mother cells and 
n the number of divisions which the multiplying differentiating 
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cells undergo, then the number of the cells capable of division at 
the steady state is 

P = mg". 
Still another possible method for maintaining a steady state is that 
P cells divide to give gP identical new cells, each of which has a 


1 1 

probability — to divide again and a probability 1 — — to differen- 
g g 

tiate and die. Equations 1 thru 5 are valid in all these cases. 


Isotope incorporation into the DNA. Two distinct cases depend- 
ing on the mode of isotope administration will be analyzed. 

Case 1. Tracer isotope is introduced into the metabolic system 
at zero time only. In this analysis we shall divide the DNA in two 
metabolic compartments: (a) Compartment A representing the DNA 
of the cells capable of division. In this compartment there is no 
cell death. (b) Compartment B representing the DNA of the matur- 
ing differentiating cells. Cell death is an important route of loss 
of isotope from the DNA of this compartment. 
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a) Kinetics of Compartment A (cells capable of division). The num- 
ber of mitoses which occur in this compartment during a time inter- 
val equal to the mean generation time is equal to the number of 
cells capable of division. If R is the total amount of DNA present 
in these cells (expressed in terms of its content of tracer element) 


: aoe 1 
and if after each mitosis the cell retains i of its original DNA and 


: 1 
acquires a fraction 1 —- 5 of new DNA, then if R, and f, are the 


amount of original and of new DNA, respectively, in these cells at 
the end of the jth generation time, we have 


Ry = R=; Ros Rens aa: 
ae | ; 
R,=R) Ri 3 Ry,=R(1—e7 7"). (6) 
i=1 


Assuming that the number of cells is very large and that there is no 
synchronization of mitoses, it can be seen that the accumulation of 
new DNA can be regarded as a continuous process resembling the 
accumulation of a stable isotope formed by a radioactive one with 
: 0.693 . ie 
a half life equal to THe I. Thus we may think of the original 
DNA as transforming exponentially into new DNA, and we may use 
the following equation to obtain RS: 


In A 


Bie Bid ao tm Yy (7) 


Notwithstanding this similarity, there is a very basic difference be- 
tween these two processes as shown in Figure 3. Figure 3a il- 
lustrates the actual accumulation of new DNA in this compartment 
in the course of five generation times, and Figure 3b illustrates 
the same process as it would have happened if it had followed the 
hypothetical model of the accumulation of a stable isotope from a 
radioactive one. Note that the total new DNA is the same in both 
cases, but that the sequence of accumulation of equal fractions of 
DNA is exactly the opposite. Thus the rate of accumulation of new 
DNA at a time 6<¢ is equal to the rate of accumulation in the 
hypothetical model at a time & = ¢- 0, and vice versa. 
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Differentiation of equation (7) gives the following expression for 


the rate of accumulation of new DNA according to the hypothetical 
model: 


b) 
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Then the actual accumulation rate at time @ is 


(ee) : es Ne ade (8) 
dt}. dt} 9 _. 6 / ; 


and the amount of new DNA Z(@) accumulated during an infinitesi- 
mal time interval d@ is 


Ink = in (t— 0) 
Se 


Z(6)=R dé. (9) 

This new DNA has the specific activity $(0) of its precursor at 
time 6; therefore, the amount of tracer isotope present in Z(6) at 
time @ is 


inf .~—28-*7,_8 
Z(0)8(0)=R —-e ! C" §(0)d 9. (10) 


If the tracer isotope is unstable and decays with a constant & 
then the amount Z*(@) which will still be present at time ¢ is 


In A 
ea 
TO" gaye kt-9) gg, (11) 


In A 
Z*(0) = Rk —e 
I 
The concentration of a tracer isotope in any given metabolic 
compartment can be represented by a polynomial function with ex- 
ponential terms (Nieset, 1955). Therefore S(6) expressed in this 
form is 


S(0) =S, e a,et, (12) 


i=l 


where “,, 7,, and v are characteristic constants of the DNA pre- 
cursor. Substituting S(6) into equation (11) we obtain 


2 In hy), 9) Y 
(0) = RS, = oe gee “yo a,c 20. (18) 
t=1 


Integration from zero to ¢ gives the total amount P* of tracer iso- 
tope present in the DNA of this compartment at time ¢, i.e., 


In A v 
Re = RS) | neler apie ha Ea ary ys 


0 i=1 
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and after integration 


Vv : es In A 
R*=RS ELS pal ES ee Gin ie (14) 
In A 


R* 
The specific activity S, of these cells is S,= Rp? which when 


substituted in (14) gives 


In A 
PMN) Be oe Ee eel | ; (15) 
é In A 


An important observation concerning equation (15) is that it is 
also applicable to the logarithmic type of cellular multiplication. 
This is confirmed by direct derivation as follows: Assuming that 
the DNA content per cell is constant and that the DNA of a cell in- 
creases / times prior to its division into A daughter cells, then if 
R.is the DNA present at zero time, the amount #&, which will be 
present at time ¢ is 

> poet 
,= Re ; (16) 

The amount of tracer element Z’(@) which is incorporated into 

the DNA at time 6 during the interval dé is 


R 1 In A 
Zotac | oo! aa eee ee eae 
dt/ 6 j : 


and the amount of tracer isotope which enters the compartment at 
that time is 
In h v 


In A 6 
Z"(6)8(8)= —— RSye 7  )> a,et dé. 


i=l 


If & is the decay constant of the tracer isotope, then the amount of 
this isotope which entered at time @ and which will still be present 
in the DNA at time ¢ is 


Z*°(0) = Z’(0)S(O)e~ #9) = 


In A In h v 
7 BSye I a,c e bhai): 


i=1 
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Integrating from zero to ¢ yields 


In 4 ae v O. : =( ka 22 8 
Teil ya ke feat | eatery 


cs ES aor as 


Substituting 7, from equation (16) into this equation and introducing 


* 
S.= = the DNA specific activity at time ¢, we obtain equation 


(15). 

In fact the previous derivation of equation (15) may be replaced 
by this one, simply by observing that the specific activity of the 
DNA of a cell capable of division is determined by the number of 
cell divisions and is independent of the ability or inability of the 
extra cells (produced by these divisions) to divide. 

Thus equation (15) can be used to determine the generation time 
I in cases where the cells multiply logarithmically, and in any case 
where the cells capable of division can be obtained devoid of ma- 
turing differentiating cells. In such cases, / can be easily deter- 


S 
mined from the value of the initial slope of the function er since 


0 
d [S, 
dt \S,/ 
z . In A\ <(4+224 
ao ee sao. as 


Be 


nh 
I 


and for ¢ = 0 


d [S, nh < 

dt (=) ! us : 
t=0 t=1 

The evaluation of J] may thus proceed as follows: The curve of 


S(@) is experimentally obtained and then graphically analyzed 
(Nieset, 1955) to obtain the constants in equation (12), which we 


S : 
may then insert in (18). The curve <a is plotted and graphically 


0 
differentiated to obtain the first derivative curve, the value of 
which at ¢ = 0 is then inserted in equation (19), which may then be 
solved for / provided that A is known (normally 4 = 2). If we wish, 
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we may take this further by inserting the calculated value of / into 
equation (15) and then apply the generalized least squares method 
(Deming, 1943) to obtain a better approximation of /. 

b) Kinetics of Compartment B (maturing differentiating cells). If 
the mean DNA content per cell is g (expressed in terms of the 
tracer element) in Compartment A, and g’ in Compartment B, then 
the rate p of new DNA formation is 


p=Aq’. (20) 
Substituting from equations (1) and (4) we have 
Mq’ D 
= os 21 
P rr L? (21) 


and 
P(g-1) , R(g-1) qv 
OS eee 


- ] - (22) 


where D is the total DNA present in Compartment B. 

Therefore the amount of tracer element 7(9) which enters the 
combined Compartments A + B of all cells at time @ during the in- 
terval dé is 

7 (0) =pd0 = SAP hl (23) 

eee. 

From this amount, an amount Z (@) will be found at a later time ¢ in 
Compartment A (vide supra). Therefore, the amount which at time 
t will be found in Compartment B will be 7(6) — Z (6), provided that 
every cell formed at time @ lives at least up to time ¢. However, 
cells do die in this compartment (vide supra); therefore, if d(¢ — 6) 
is the probability of survival of these cells up to time ¢, then the 
amount 6(6) which will be found at that time is 


5(8) = [7(9) - Z(A] d(t - 4), 
and the amount of tracer isotope will be 
5 (8)S(6) = [7 (0) - 2 (A)] d (¢ - 6)S(6). 


However, if the isotope decays with a constant &, then the amount 
of tracer isotope 5*(6) still present at time ¢ is 


5*(9) = [7 (0) — Z(A)] b(t - 0)8 (A)e~*Ft-®, 
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Substituting 7(6), Z(6), and S(6) from equations (23), (9), and (12), 
we obtain 


e 


R a inal =¥ 
SAE) ar 28s le-9 £- inne i= J ee ee 


Vv 
y. ae! p(t - ade. 


integrating from zero to ¢, we obtain the total amount D* of tracer 
isotope present in the DNA of this compartment at time ¢ 


R Vv t Ge wt 8) 
t= — § a, f jig - 1) © in he |. 
/ oD at A q 


Bhp Ft) 


b(t - 0)d6. 
Placing ¢ - 6 = z when d6@ = —dz, we obtain 


uv q’- aR 
tao hes: ee f @-) £- inne 4 iE 


(24) 
pakke T dx d (x)dz. 
From equations (21) and (22) we derive 
i q 
EAA ATC AY a 
D* 
Substituting in (24) and placing Sai S, we have 
v t 2. _inh : 
Pia tS ey ae" [ l-» ea Ree i x 
8) L tee 1)q iz fy) q (25) 
nae NE d (x)da. 
Differentiating in respect to ¢, we obtain 
v , _inh 
di)! tI, pot mis ]ee. 
t 
0 (g-He q (26) 


In h 


L f. y — x 

>: are" f j- 9 L-in ne i ‘le baal sto} 
q 

i=1 . 
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This equation may now be rearranged to give the value of one of 
the integrals in terms of the others, i.e., 


oink 


t Ke aaa EN = r.)x 
aye f [i= 1) © - in he £ fas O* b (x)da = 
q 
0 


—1\ i el & v , ao ale 
BS) PA Eada ONENESS o io-) £- inne ml 
fs q dt \S, q 


t=1 


k Vv t gy’ -Ay 
—kt rot 
€ Os arse | li@=1) £ mae |. 


0 


asthe (oe| : 
Substituting this into equation (25), we have 


ne Geuk. He | OS ee ae) 
ee (9-1) 7 — In he x (27) 


0 


—(k+r,)x “0; q” Sie a ee 
e i $ (a)de - (g—1)—-Inhe .* Je ater 
re: q 
Repeating this differentiation and back substitution vy — 1 times 


more, we can eliminate all the integrals and obtain a differential 
equation of the form 


S, v d‘ 8, 1 dg (2) 
Serle L (g-1)¢ a Dail Saree camer os! 


where c,, C,,-..,¢,, are known constants, and fis fo(s+++ sf (4) 
are known minctond at t. 


Thus if vy = 1, it can be seen from (27) that 


1 x , ee at 
oy=-— and f,()=-—! \(g-1) H-inhe 7 ‘| em, 
r r q 


1 1 
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hr Of, 


2 , inh, 
no--{[p a ies “a)| fo- 2 -ina : |. 
i= 1 2 y 


1 1 1 
When v = 2, thence, =-|— + —], ce, = : 


and 


2 aX q’ ae: a ¢ 
fa(t) = )° Pn te : —nve: ew kt, 
172 


i=1 


Similarly if vy = 3, we obtain 


1 1 ih ih i! 1 
fe argue ts Pa foe fa fs if o Fs 
1 = 
oF age ’ f,()=-(,+¥,+V,)e die 
Pyrat; 


f,(2) = (wy, 7 2u,)e** 


8 , Ink 
Oo; vl Te-] i —kt 
ek ) -1)—-lInhkhe e 
f,(é : rr? ( ) q ’ 


3 


and 
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, _inA 
lo-» inne I ak 


: a; In A -A, 
Hie Sipe pees f-na(ee 4). ° 


The boundary conditions of equation (28) are ¢=0, d(t) = 1; t=~, 
¢(t)=0. Now it is reasonable to assume that a minimal but not 
infinitesimal length of time elapses from the formation of a cell to 
its autolysis. During this interval there will be no cell death (ac- 
cording to our definition of L), i.e., d(¢) will remain constant and 
equal to unity and all of its derivatives will be zero. Thus, the 
initial conditions at ¢ = 0 are 


Tis se] Ee) P peed awit 


t=0 


and 


Introducing these values in equation (28), we obtain for ¢ = 0 
Siege fe) Senn f (0 (29) 
Ct Sel es 29 
= * dt*\S, 7 a (9g - AF 


Differentiating (28) in respect to ¢ and taking the value of its de- 
rivative at zero time, we obtain 


d/s vil 
dt ie) “ ‘2 Ci ie fae) = a pee Ls ) ss (30) 
alice <4 ati\So/,, Lg-el a |), 


1=2 
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Dividing (29) by (30) we have 


v t 
ace Sa 
att \s 
t=0 


i=l 0 


eet On 


d /s vt af] ’ 
2 See 
dt \s, i-l Geils dt 

O/ t20 


RaQ. f tnd JAM 


The left-hand member of this equation can be obtained as follows: 
Following the derivation of equation (12), we derive equation (28) 
and express the constants c, in terms of the constants r;. The 
functions f;(¢) are also derived and expressed in terms of ti Oy 7, 


S 
and t. The derivatives of —* can be obtained by successive graph- 
0 


: ore Ss 
ical differentiations of the experimentally obtained curve —¢. Fi- 


nally, taking the values of these derivatives at ¢ = 0, finding the 
values of f,(¢) and d/,(t)/dt for ¢ = 0 (provided that the constants 


A, g, and Be are known), substituting into (31), and solving, we ob- 
qg 


tain 7. L is then found by substitution into (29) or (30). Knowing 
L and /, equation (28) can be solved by the usual methods (Margenau 
and Murphy, 1956) to yield the probability of cell survival ¢(¢). 
The probability of cell death prior to a time ¢ is 1- ¢(¢), and the 
probability density of cell death, or the rate of cell death, is 

dd (t) 


d 
q -eOl-- FF bles 


Thus graphically differentiating the curve ¢(¢) and plotting the 
_ dd (2) 
dt 


, we obtain the curve of the probability of cell death, or of 


the rate of cell death. The median life span L is the time when 


is maximum, 


~ d 
¢ (t) = 0.5, and the mode life span L is that when — 


c) Kinetics of the Combined Compartments A+B (all cells). If Q 
is the total DNA of all cells, and Q* the isotope tracer present in 
it, then 

Q=R+D and Q*=R*+D*. 
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Taking R* and D* from (14) and (24), we have 


ek ot eae In h | gear) 
Q*= 7S, )) ae Ta ee + 


(32) 


, 


: q -inh x =—(khit rs) x 
[ eae Peet i é i” $ (x)dz 
0 


R 
Combining equations (5), (20), and (22) and solving for 7? we obtain 


| x 


Tq 


2a S 33 
I (g-1)Le+l (33) 
However, by definition, T = P+ M. Therefore, Q = Pg+Mq’. 
Substituting M from (2), P from (3), and solving for T, we have 
os Ql(g-1)L+7) 
ql+(g-1)q’°L° 
Placing T into (38), we obtain 
R Q 
—= Std pate ae (34) 
I (g-lyr’le+d 


R 
Finally substituting T in (32) and representing the specific ac- 


ae he 
tivity oy. by So we have 


So. esac see Ini 
S, (g-l¢Lsal ae Ink 
k+?,; + — 
I 
Ink t y ink. 
= meal lo) © ae |x @s) 
0 


e Rtr is 


p (z)dxr , 
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or 
Sq q 
— + F(t) = ——_—_____ x 
S -1¢7L i 
0 (g-1)¢’L+q (36) 
: Sis is “a -in4, —(k + 
De [ (g-1)—-Inhe !? |- 'O* % (@)de, 
i=1 0 q 
where 
ae | one e r+ inh 
F (t) = oe eT (k+r,+ ae 
G-De+ go ink 
I 


Following the same procedure which we used to derive equation 
(28) from (25), we can derive from (36) a differential equation of 
the form 


S, : d‘ /S,\ d'‘F(d) 
gar O*). « [Fa(34)+ wl - 


a (37) 
q a *6 (t) 
ee Oar dito 


where the constants c, and the functions /,(¢) are identical to those 
of equation (28). The boundary and initial conditions are also the 


S 
same as those of equation (28). Therefore since a = 0 and F(t) =0 
0 
at ¢ = 0, equation (37) at zero time becomes 


y S d'‘F (t) q 
2 es E (=) ME EF. | Ax D (g-1)¢Liq 1,(0)s 


i=1 
te eT 
(g-1)¢L+qQ 


and by placing F (¢) = u(t), we have 


Vv d! S 
ya ast) 
7 a O/ t=0 


(38) 


q _< oe }. 
(gil) 9 D+ 5 {ho dL “ dt? J. .4 
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Differentiating equation (37) and taking the value of the derivative 
at ¢ = 0, we obtain 


d s v git! S, 
ae) + Do cgi ise ti 
0/t=0 1 0/ t=0 


i= 


q site] Ber. a } eu. 
(g-1)qLiql [e ; a i dtt*} eae 


Dividing (38) by (39), we have 


(40) 


: dé u(t) 
f,(0)- )) c.| =e | 


i=1 =" 
v git} t : 
Ee @) = SE c, = I 
dt scneGee at enc 


The left-hand member of this equation is evaluated the same way 
as that of equation (31). The right-hand member is a function of 
the generation time I and can be used for its determination. L can 
be found by substituting / into equation (38), provided that the con- 


stants h, g, and 5 are known. Equation (37) now can be solved to 
q 


give the function #(¢), from which by graphical differentiation we 
dd(t 
can obtain the probability of cell death -— 8 : 


Case 2: The specific activity of the DNA precursor is maintained 
constant. This is a special case of the more general Case 1, for 
which we have S(6) = S, = S = constant. Therefore the equations 
which are applicable in this case can be obtained directly from 
those of Case 1 by placing v= 1, a, =1, andr, =0. 

a) Kinetics of Compartment A, and of logarithmic multiplication. 
When the above-mentioned substitutions are made in equations (15), 
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I 
(18), and (19), we obtain 


Inh 
S, Inh E 3 oo Met 24 


Ss ¥ kl+ilnhk (41) 
d [S; Inh -(x+424), 
+ (z)-"" 4 eee 
and 
d [{S, Inh (43 
Ss Re ay 


Thus, following graphical differentiation of the experimentally 
S 
obtained curve 3? equation (43) is used to determine 7. The gen- 


eralized least squares method may then be applied to yield the best 
value of /. 

b) Kinetics of Compartment B. The same substitutions when made 
in equations (25) and (26) yield: 

By 1 q 


t q’ —inh, k 
Te oor jo-92 - inne |: p(x)dx,. (44) 


and 


_inh, 
£ (%) ' = be fae : Jer 4, (45) 


At t = 0, d(¢) =1, then 


d (*) tae, E ee (46) 


a\s},., EL. @-a 


Assuming that the constants h, g, and 2 are known, the mean life 

span can be calculated by entering these constants and the zero 
: S 

time value of the graphically obtained first derivative of e d 


Differentiation of (45) leads to 


#(%)- + {{s- q inh oes pea) ie 
d?\s/ L (g -1)¢ at (47) 


nh 
Co ard are aa aa 
i= 
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_ aot) ; 
Att =0, d(¢) = 1 and taking ie 0 (vide supra), we have 
2 In h Inh 
ee 1 ania eer mi 
di? \8 Jo 58 ub hg ae I 
Dividing (46) by (48) and solving for /, we obtain 


d [84 
q (In h)? — (= 
dt \S/ 20 


[(g - 1) 9’ - gin’ s (=) mo Soy, ue 


The mean generation time | thus can be obtained from equation (49) 
in which the zero time values of the first and second derivative 
curves, obtained by two successive graphical differentiations of 


i= (49) 


S 
the curve re are used. 


d {8S 
Knowing L, J, and — “) , we can calculate the probability of 


dt \S 
cell survival from (45) which gives 
Le*t @ ‘a 
dt \S& 
$ (2) = fei Se ee ea ye (50) 
“on AY. saat 
(9 -1)q’ 


: d 
We may then obtain the probability of cell death — a8. From 


the distribution of deaths, the median, and mode life spans can be 
found. 


NY 
In regard to the shape of the curve - , we observe from equation 


, d 
(45) that, since />0 and 0< g(t) <1, ai (3) is always greater 


than zero (¢ > 0) and has an asymptote parallel to the ¢ axis. Thus 
S 
the curve “e does not have any maxima or minima but increases 


steadily possessing an upper asymptote parallel to the ¢ axis. In 

d¢(t) 
dt 

that some of the terms in equation (47) are always positive and the 


addition, because 0 < J(t) <1, 


£0, and g >1, we observe 
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: d2a(/S 
rest are always negative. Thus 770) (=) can be either positive or 


negative depending on the function ¢(¢), and equation (47) may 
a 8 
have one or more positive roots. Therefore the curve 5 may pos- 


sess one or more inflection points (see for example Figure 6). 
c) Kinetics of the combined Compartments A +B. When the proper 
substitutions are made in equation (35), we obtain 


a ee ee 
S (g-l)q’Ligql Inh 
k+ — 


: (51) 
t , -inh, 
} lo-9 © - mre I et o(e) de 
0 q ; 


A special form of this equation is the simple equation which we 
have previously derived and used (Osgood, Seaman, Tivey, and 
Rigas, 1954). This equation can be obtained from (51), assuming 


Us 


(2) = 1 (i.e., no cell death), ‘A = 2, g =2, and =e 1. Integration 
q 


then gives oie ae 
Sebtedddies De: 
However the use of this equation is very limited because of the un- 
realistic assumption that d(z) =1. 

Two successive differentiations of (51) yield 


Pai a 
dt \S (g-1)q°Li+ql 
, _ink - (n+ 424): (52) 
{[o-1 © - inne : ‘Jem oto inne! : \, 
q 


d* (Sq\ _ q ie 
dt? \s}/ (g-lq’Lt+ql 


, Ink; 
{[@-v © = mae? eee (53) 


ich (.- ae ah A oH g(o}. 
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At t = 0, 4(t) = 1; therefore from (52) we obtain 


@ 7) rn aes Se (54) 
dt\S/,.5 g-DaL+ql 


dt? \S 
However, we observe that since lim d(¢) = 0, then for ¢ > r equa- 
t>T 


d2 
In this case it is not possible to obtain / from —> (5) : 
t=0 


tion (52) becomes 


d (= qlnh _ ts 
t>T 


dt \S. ofg = De Lgl 


Substituting the expression (g ~ 1)q¢’L + q/ from (54) and solving 
for | we obtain 


In h 
eee Baek a (55) 


d [S ; 
@-o9 § (2) 
poe ay ie See RT ee 
ad {8S 
ma 
dt MS ies 


where ¢>~+7. We are led thus to consider the function 


In 


Inhkt 


I(t) = - (56) 


for which lim /(t) =/. 
t>T 


It can be shown theoretically that this function is zero at zero 
time; it has an asymptote parallel to the /(¢) axis at a time 6(0<6<7); 
and it is negative for ¢ < @ and positive for ¢>@. Thus the mean 


, 


generation time I can be obtained, if A, g, and Ee are known, by 
q 


calculating and plotting the values of this function as shown in 
Figure 5. Another method which can be used to derive the values 
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of L and /, if the percentage of cells capable of division p = 100 a 
At 


and the constants g, and :® are known, is to solve equations (3) 
q 


and (54) as a simultaneous system. Equation (52) may now be 
solved to obtain 


d [8S —inh, 
l(g -1)q’L+qle* = me -~qlnhe ! 
$(¢) = ————___*_______._ 657) 
(g-1)q’~qln yates 
d $(t) 
dt 
mode life spans, are obtained as in the previous sections. 


The probability of cell death ~ , and hence the median and 


S 
In regard to the shape of curve ee we observe that since / > 0, 


and a(t) 
dt 


g>1,0< ¢(d <1, <0, equations (52) and (53) have no 


S 
roots other than ¢ = ~, and therefore the curve = does not possess 


any maxima, minima, or inflection points, but continuously in- 
creases, having an asymptote parallel to the ¢ axis as shown in 
Figure 4. 


Nonrandom miztures of cells. It is conceivable that sampling 
cells from tissues may yield a mixture of cells, some dividing and 
some maturing differentiating, in which the ratio of cells capable 
of division to that of maturing differentiating cells may not be 
identical to the ratio in the tissue. Thus if f is the ¢rue ratio in 
the tissue and y is the apparent ratio in the sample, then the amount 
E of DNA (in terms of tracer element) of a sample containing all 
the maturing differentiating cells would be 


and its tracer isotope 


| 
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Substituting from (14) and (24) we have 


lnk 
- ryt mE: | 
B* = — Inh ie ne 
So i at aoa 


Ble haar a 


te matt gs S “ix 
> oe jo-»" —iInkene | a "iin d(x) dz 
q 


i=1 0 


R 
Equation (34) gives Fi in terms of Q, which in turn can be expressed 


in terms of E as follows: Since Q=R+D and D=E-yR/f, it 
follows that 


bs El(g-1)q’Li gl] 
, bd 

-1 L+=ql 

(g-l)¢ bmi 


Substituting in equation (34), we obtain 


R Eq 


ae 
(g-1)q°L+—ql 
B 


KE * 
Substituting this in the equation for E* and placing — =S,. we 


have 


—~ 


e g 


—_- =: _—_—_—_—_— —_——_ XxX 
S 


to oly 
oe Cela 


7: 
Inh 
v ~ me Ink 
tee a Guan jr-« Csreted rt (58) 


pat kre 
I 


. q’ ah ae 
if tg tml) alse I e "i * b (ax) da 
ils 


It can be seen that this equation is the generalized form expressing 
the kinetics in the steady state, since when y =f it becomes 


identical to (85), and when y = 0 (cells capable of division being 
absent) it becomes identical to (25). 
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Following the same procedure that was used in the derivation of 
equation (37), we obtain from (58) 


S : fe S ago 
—£ + G(t .|—{(— = 2 
pa eke) de “tL ai Fes) * at! 


0 i= 


Ho (59) 
q d'~' df(t) 
y 2, file a eas 
(g-1)q°L + —q/ im! 
B 
where 
ad rt 
" taeda e+ oa 
SUG Fa eta eae 
ae een eee a ie | 


ba ; 
(g-1)¢°E+ — ql i*1 ker, + — 
g-1)49 a4 aes 


c, and f,;(¢) are the same as in (37). 

The special equations in cases in which the specific activity of 
the precursor is constant can be easily derived by proper substitu- 
tions as in the preceding section. 


Illustrative example. As an illustration of these methods, we 
shall analyze the data which we have obtained on the in vivo up- 
take of P®? into the DNA of the blood leukocytes of a patient with 
chronic leukemic granulocytic leukemia. This patient’s blood leu- 
kocyte count remained constant within experimental error at the 
level of 315,000 WBC per cmm. The differential counts indicated 
that 95 per cent of the leukocytes belonged to the granulocytic 
series, the rest being mainly lymphocytes, and that approximately 
20 per cent of the granulocytes possessed nucleoli. The leukocytic 
tissue of this patient was thus considered to fulfill the requirement 
of the steady state, and therefore to be suitable for this study. The 
plasma inorganic phosphorus was assumed to be the DNA precursor 
and its specific activity was maintained constant by calculated 
daily injections of P?? (Osgood and Tivey, 1950) following an 
initial intravenous injection of 200 pe of P®? in the form of sodium 
phosphate. Leukocytes from blood samples obtained periodically 
thereafter were separated by the phytohemagglutinin technique (Li 
and Osgood, 1949}; Rigas and Osgood, 1955) and their DNA isolated 
as previously described (Rigas, Duerst, Jump, and Osgood, 1956). 
Subsequently, the specific activity of the DNA phosphorus was 
determined and was expressed as a fraction of the constant specific 


wa 
ra 
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S 
activity of the plasma phosphorus, i.e., as rn These and other 


pertinent data are tabulated in Table I. 

In Figure 4 we have plotted the data of the specific activity of 
the DNA phosphorus against time, and a smooth curve (a) is drawn 
through them. The next step involves the graphical differentiation 
of the curve (a) which yields curve (b). The apparent absence of 
inflection points from curve (a) suggests that our leukocyte samples 
contained cells capable of division as well as maturing differenti- 
ating cells. This conclusion is supported by cytological evidence 
indicating that leukocytes of the granulocytic series which possess 
nucleoli are capable of division (Osgood and Brooke, 1955), and 
such cells were present in these samples to the extent of 20 per 
cent. Thus, assuming that the ratio of the cells capable of divi- 
sion to the maturing differentiating cells as found in the blood sam- 
ples is equal to the true ratio of these cells in the body of the 
patient, i.e., B = y, we can analyze the data according to the ki- 
netics of the combined Compartments A+B. Taking the zero in- 


, 


tercept of curve (b) and assuming that g = 2 and an 1, we obtain 
q 
from equation (54) 


1 1 
L+]= ——— = — or L +] =4.3 days. 


d (S, 0.23 
dt \S/,., 


Equation (56) may be used here to calculate / (using the decay 
constant of P°?, % = 0.0485). Values of /(¢) have thus been ob- 


A 


; ; d [S 
tained, taking the values of ai (*) from curve (b). These values 


TABLE I 

: Cells of Leukocytes ; 
Time Leukocytes granulocytic with Specific KY 

in per cmm of pF a wunlesli activity 
days blood x x of plasma 
0 306,750 96.8 14.0 - 0.00 
0.8 306 ,000 96.4 15.0 0.030 0.19 
2.8 302,000 91.2 25.4 , - 0,027 0.58 
4.8 312,000 95.8 24.2 0.027 0.74 
6.8 350,000 93.0 20.8 0.027 0.78 


Mean 315,350 94.6 19.9 0.028 
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<3 (0); OR a (22). OR 


(d). 


dA(1) 
dtr 


Pir), (c) i OR - 


FIGURE 4 


are tabulated in Table II and have been plotted in Figure 5. The 
limiting value of /(¢), i.e., 7, cannot be placed with much accuracy 
since data were not obtained for a sufficient time for the curve of 
I(t) to level off. Extrapolation of /(¢) suggests that / is about 1 
day. However, we can apply the alternative method of calculating 
I from equation (3), assuming that the cells with nucleoli are the 
cells capable of division. Since L +/]=4.3 days and the mean 


= 
percentage of these cells is p = 100 mp 20, equation (3) gives / = 


0.9 days which is essentially the same as that indicated from /(t¢). 
However, because of the relatively large errors involved in doing 
differential leukocyte counts, this figure may not be significantly 
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TABLE 

Time d [Sq =a do (bt) 

in days dt (52) A, 9 (2) dt 
0.5 0.227 — 0.946 0.998 ~ 
1.0 0.220 - 1.931 0.989 0.030 
1.5 0.212 — 3.007 0.972 0.076 
2.0 0.196 — 4.784 0.912 0.145 
2.5 0.170 — 9.967 0.800 0.242 
3.0 0.137 114.231 0.645 0.400 
3.5 0.090 5.866 0.424 0.326 
4.0 0.060 3.960 0.320 0.230 
4.5 0.045 2.946 0.215 0.173 
5.0 0.029 2.353 0.142 0.1384 
5.5 0.017 1.918 0.082 0.090 
6.0 0.011 1.734 0.052 0.060 
6.5 0.007 1.591 0.034 0.040 


different from one day, and for simplicity we shall assume 7 = 1.0 
day. The mean life span of the leukocytes of this patient would 
thus be L = 4.3 — 1.0 = 3.3 days. 

Values of the probability of cell survival d(¢) can be calculated 
with the aid of equation (57). These values have been tabulated 
in Table II and plotted in Figure 4, curve (c). From this curve it 
can be seen that there is practically no cell death during the first day 
and that 50 per cent of the cells are dead at 3.35 days. This point 


FIGURE 5 
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is the median life span L of the leukocytes. Graphical differentia- 
d $(t) 

dt 
(d) of Figure 4 which gives the probability or rate of cell death at 
different cell ages. The maximum ordinate of this curve gives the 
mode life span L which in this case is 3.2 days. 

In order to demonstrate the analysis of a case identical to the 
preceding one with the exception that the leukocyte samples would 
contain only maturing differentiating cells (when the kinetics of 


tion of d(¢), i.e., of curve (c), and plotting of - yields curve 


S 
Compartment B would apply), the theoretical values of Fe have 


been calculated, assuming L =3.3 days, ]=1.0 day, and ¢(t) is 
identical to curve (c) of Figure 4. These values have been tabu- 
lated in Table III and plotted in Figure 6, curve (a). This curve, in 
contrast to that of Figure 4, has an inflection point at about 1.75 
days. Starting from this curve, we shall derive the same figures as 
in the previous case. 

Graphical differentiation of this curve yields curve (b) of Figure 
6, i.e., the curve for - * . Curve (b) is again differentiated 


S 
Figure 7. Taking the zero intercepts of these last two curves and 
applying equation (49), we have 
0.693? x 0.093 
~ 0.307 (0.14 + 0.0485 + 0.093) 


; _ a (8, 
graphically to yield the curve wa \s)” which has been plotted in 


I = 1.006, 
and rounding off, / = 1.0 day. 
The mean life span L calculated from equation (46) is 


0.307 


Ye: pitguiwa aes 
0.093 nh 


Values of ¢(t) can be calculated now with the aid of equation (50). 
The resulting values are shown in the last column of Table III and 
have been plotted in Figure 6, curve (c). Curve (d) of this figure 
can be obtained in the same way as that of Figure 4. 

Remarks. The preceding mathematical analysis illustrates the 


4 


necessity of knowing the value of the constants h, g, and &: as 
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TABLE OI 

i Ss a\s, a2 

in days S at \8S dt S 
0.5 0.063 0.152 0.089 1.000 
1.0 0.148 0.186 0.051 0.987 
1.5 0.247 0.207 0.019 0.970 
2.0 0.352 0.207 —0.017 0.913 
2.5 0.451 0.189 —0.053 0.801 
3.0 0.538 0.154 —0.084 0.643 
3.5 0.600 0.102 —0.080 0.426 
4.0 0.644 0.076 —0.055 0.319 
4.5 0.675 0.051 — 0.038 0.216 
5.0 0.697 0.033 —0.027 0.142 
5.5 0.711 0.019 —0.021 0.084 
6.0 0.720 0.012 —0.015 0.054 
6.5 0.725 0.007 —0.010 0.032 


well as whether the cells present in the samples are cells capable 
of division (Compartment A), or maturing differentiating cells (Com- 
partment B), or a random or nonrandom mixture of the two, so that 
the appropriate equations can be applied in the derivation of the 
life span and/or the generation time of the cells. This information 
generally should be obtained from cytological and histological 
studies of the tissue under consideration and of the analyzed cells. 
However, knowledge of the absolute values of g and q’; i.e., of the 
mean DNA content of each of the two types of cells, is notnecessary. 
Furthermore, the exact ratio of the two types of cells present in 
the samples should be known only in cases where nonrandom mix- 
tures of cells are analyzed. Otherwise it is sufficient to know 
whether the samples are random mixtures or devoid of one or the 
other type of cell. When the samples are random mixtures or devoid 
of cells capable of division, one can calculate both L and J from 
the DNA specific activity curve and then determine the proportion 


FE 
of these cells a utilizing equation (3). Only the mean generation 


time I can be calculated, though, if the samples consist entirely of 
cells capable of division. It should be emphasized again that as 
far as this theoretical analysis is concerned, it is quite immaterial 
whether the cells capable of division are identical to each other as 
in Figure 1 or fall into two categories as in Figure 2. Furthermore, 
this analysis is applicable whether each dividing cell always gives 
two daughter cells, one of which is capable of division, while the 
other is a maturing differentiating cell (Osgood, 1953; 1954; 1957), 
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FIGURE 6 


or whether the total number of cells capable of division yield g 
daughter cells, each of which has equal chances to either divide 


1 
again or differentiate and die, with the result that only — of them 
g 


(regardless of origin) divide again to maintain the steady state. 
Thus the equations derived here are valid no matter which one of 
these mechanisms is responsible for the maintenance of the steady 
state. 

An important point in this analysis is that if the steady state is 
maintained through a mechanism similar to that shown in Figure 2, 
then the life span L, as defined in this paper, does not include the 
interval from the time of the formation of a multiplying differenti- 
ating cell to that of its last division. Thus the magnitude of this 
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FIGURE 7 


interval cannot be obtained from data on the isotope incorporation 
into the DNA unless additional information is available in regard 
to the number of cell divisions n which a multiplying differentiating 
cell can undergo before its progeny is converted to maturing dif- 
ferentiating cells. 

The plausible probabilistic features of 7, A, and possibly of n, 
should not be overlooked. However in the steady state, there is a 
constant rate of cell division and of new DNA formation and flow 
from Compartment A to Compartment B, due to the balance of: (a) 
the cells which divide ahead of time by those which have delayed 


1 
division, (b) the ones which retain more than + of their original 


DNA by those which retain less, and (c) the ones which miss a 
number of divisions by those which divide an extra number of times. 
Thus, although it is possible to determine the mean generation time 
I from the rate of isotope incorporation into the DNA, it is impos- 
sible to gain any information on the distribution function of / in re- 
gard to individual cells when the pooled DNA of populations of 
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cells is studied. This is also true in the case where the generation 
time changes in an orderly fashion as the differentiation of the 
multiplying differentiating cells proceeds, as has been suggested 
for the hemopoietic tissue (Patt, 1957). Here again, as long as a 
steady state is maintained, the specific activity curve of the pooled 
DNA is not affected by the variance of / among individual cells. 
Thus, the generation time that is derived from such studies is a 
mean figure, ‘the variance of which and its possible correlation to 
the stage of cell differentiation can only be obtained from cytologic 
or other studies of the individual cells. 

The steady state requirement for the applicability of this method 
is an important one. Although an absolute steady state may never 
exist in the tissues of an animal, we may assume the steady state 
as being a reasonable approximation whenever the rate of growth 
or regression of the tissue is negligible in relation to the rates of 
cell division and cell death. As mentioned previously, this analysis 
necessitates that the steady state be maintained by cell division 
and cell death, the rates of which are equal and constant. It is 
conceivable, ‘however, that a steady state can be maintained also 
if these two rates change with time but remain always equal to 
each other. In this kind of steady state it can be seen that since 
the cells which are formed at any instant die on the average L time 
later, the functions representing the changing rates of cell division 
and cell death are out of phase in respect to each other by the time 
interval L. To balance each other, they should therefore be identi- 
cal periodic functions of time, while the lag of phase of the death 
rate should be an integer multiple of their period and should have 
a negligible variance. Thus a steady state that would be main- 
tained by such a mechanism appears to be a unique case and, 
therefore, rather improbable. 

The main advantage of this method of analysis over others previ- 
ously proposed (Kline and Clifton, 1952a,b; Ottesen, 1954; Osgood, 
Seaman, Tivey, and Rigas, 1954) is that it provides an exact method 
of analysis of data even in cases where the rate of death of newly 
formed cells may be quite large following the time of their forma- 


S 
tion, i.e., whenever the curve —* does not show the necessary 
) 
plateau to separate the function ¢(¢). It also enables the determi- 
tion not only of the mean, mode, ‘median, and distribution of the 


life span of the cells, but also of their mean generation time. 
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A somewhat different approach to the principle of biotopological map- 
ping, discussed in previous publications, is given. The organism is con- 
sidered as a set of properties, each of which is in its turn a set of numer- 
ous subproperties which are logically included in the corresponding 
properties. Topology is introduced by an appropriate definition of neigh- 
borhoods, and four postulates are stated which concern the mapping of 
the spaces corresponding to higher organisms on those of lower ones. A 
number of conclusions are drawn from the postulates. Some of them cor- 
respond to well-known facts. For example, in man and some higher organ- 
isms appropriate emotional stimuli should produce gastrointestinal or 
cardiovascular disturbances; or some microorganisms should produce sub- 
stances harmful to other microorganisms (antibiotics), Some other con- 
clusions are still awaiting verification. One of them is, for example, that 
there must exist unicellular organisms which produce antibodies to ap- 
propriate antigens. 


In a series of publications (Rashevsky, 1954, hereinafter re- 
ferred to as I; 1955a—II; 1955b—III; 1955c—IV; 1955d—V; 1956a 
—VI; 1956b—VII; 1956c—VIII; 1957a—IX; 1957b—X), we have 
suggested and outlined a topological approach to general biology. 
The aim is to emphasize the relational aspects of biology and to ex- 
press the unity of the organic world. This approach is in no way in- 
tended as a substitute for the biophysical, metric approach used 
hitherto in mathematical biology. The two should rather supplement 
each other. Great as are the advances of the current metric physico- 
mathematical approach to biology, they leave entirely outside of 
their scope a large number of important relational phenomena. They 
also do not express the basic similarities between the most dif- 
ferent organisms. Such phenomena as sensitivity to stimuli, some 


72 N. RASHEVSKY 


forms of motion, digestion, excretion of wastes, etc., are present in 
one form or another in every organism, unicellular or multicellular. 
The quantitative aspects and frequently the physicochemical mech- 
anisms of any of these properties are very different in different or- 
ganisms. Yet they not only are present in all organisms, but the 
basic relations between those properties are the same for all 
organisms. 

This general observation leads to the formulation of the principle 
of biotopological mapping, as discussed in I. In some way all or- 
ganisms can be mapped on each other and all of them on some con- 
ceivable simplest organism, which may either exist or have existed 
or which may be conceived as an abstraction possessing the mini- 
mum number of properties which makes us recognize an organism as 
such. For lack of a better name, we designated this simplest or- 
ganism as the primordial. Ifthe relational structure of the primordial 
organism were known and if the mathematical nature of the mapping 
were known, we would be in a position, barring purely mathemati- 
cal difficulties, to describe the relational aspects of all organisms. 


The purpose of the present paper is to give a somewhat more 
general and yet a somewhat more precise formulation of the princi- 
ple of biotopological mapping and to show how a number of known 
biological facts follows rather directly from this formulation and 
new phenomena are predicted by it. 


In I we pointed out that an organism is a set of what we called 
‘*biological functions,’? connected with each other by certain rela- 
tions which impart on that set an ‘‘organization.’’? Examples of 
‘*biological functions’’ are sensitivity, locomotion, digestion, se- 
cretion, etc. We say, for example, that the function of the stomach 
is to digest, the function of a gland is to secrete. In a mathemati- 
cal discussion we have to use frequently the word ‘‘function’’ in 
an entirely different sense, namely, as the mathematician under- 
stands it, e.g., x= f(y). Though in I we agreed to use the word 
function in the mathematical sense, while explicitly using the com- 
bination of words ‘‘biological function’’ in the sense in which the 
biologist uses it, yet the notation is not a convenient one. Some 
‘‘biological functions’’ are characterized by their intensity, for ex- 
ample, the intensity or rate of secretion. Situations may arise 
where we would have such monstrous expressions as ‘‘the inten- 
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sity of a biological function is a function of the intensity of another 
biological function.”’ 

We shall therefore henceforth change our previous terminology and 
denote such things as sensitivity, locomotion, digestion, absorp- 
tion, secretion, etc., by the word “‘property’’ of an organism. We 
shall denote different properties by the letter P provided with a 
capital Latin subscript to designate the specific property. Thus 
P, stands for sensitivity to external stimuli; P.—for property to 
conduct excitation; P,—for orderly movement (as contrasted with 
random thermal agitation); P,—for absorption, etc. Lower case 
Latin letters will be used for running subscripts. Thus we shall 
write P,, (i=S, C, M...). 

In a unicellular organism the one cell, which constitutes it, pos- 
sesses all the basic properties P;. Those properties are related in 
a manner to be discussed below. One of the characteristics of 
multicellular organisms is that different cells specialize in dif- 
ferent properties. Some cells retain some properties and lose 
others. The relation between two different properties P; and P,, 
whether they belong to the same cell of a multicellular organism or 
to different cells, remains the same as in a unicellular organism. 
For example, one type of cell may be sensitive to external stimuli, 
but may not be motile; another cell may not be sensitive but may 
possess motility. Yet just as in a unicellular organism a stimulus 
under certain conditions produces a movement, so will a stimulus 
applied to the sensory cell result in a movement of the motile cell. 

Another characteristic of more complex multicellular organisms 
is the much larger number of different properties as compared with 
the simpler or unicellular ones. Thus a sensory cell in a lower 
organism may be sensitive to chemical stimuli in general. In a 
higher organism, different sensory cells may be sensitive to dif- 
ferent chemical stimuli. For example, different taste cells of the 
human tongue are sensitive to salt, sour, sweet, and bitter. In I 
we designated those additional properties as ‘‘subsidiary biological 
functions.’’ In our present terminology it is more appropriate to 
speak of ‘“subproperties.’’ The sensitivity to bitterness, which we 
shall denote by Ps, and the sensitivity to sourness P,, are in- 
cluded in the sensitivity P;. Thus P, is a set of subproperties 
Peano me lice «le, Be Peg 3 are subsets of P,, or elements of P,. 
The same holds for any other property P,. Denoting a ‘‘subprop- 
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erty’ by P,, (i =S,C,M...; 4=1,2,3.. .)(in other words numbering 
the subproperties in an arbitrary manner) we have 


EE ng (1) 


When we examine the mapping of different organisms on each 
other, we see that the set S,(P,,) of the subproperties of P; in one 
organism maps on the set S, (Pig) of another organism; and both 
map on the set P, of the primordial. The set P; is thus the largest 
or the most inclusive of all subsets P, ,. 

The sets P, may be conceivably infinite. Thus P; contains amongst 
its subsets the sensitivity P,, to radiation of wave length A. Thus 
P, has, at least theoretically, the power of the continuum. Actu- 
ally, however, in any given concrete organism, the sets P, are all 
finite, though some of them may have very large cardinal numbers. 
An arbitrary ordering of the subsets P,,, as done above, does not 
present therefore any practical difficulty. 

It must be also noted that two subsets P,,, and Pig of P; arenot 
necessarily disjoined, so that we may have 


Pig 0 Pig = 0.A.P.40 Pg #0, (2) 


where 4 denotes the exclusive ‘‘or.’” As an example we have the 


secretion of proteins and the secretion of regulatory substances. 
Both are included in ‘‘secretion’’ in general. Some secreted pro- 
teins are regulatory substances, and some regulatory substances 
are proteins. 

Let us examine more closely some of the sets P,. 

We have already discussed some of the subsets P;, of P;. Asa 
further example let us consider that a painful stimulus may be a 
stimulus of any modality with an intensity above a certain in- 
jurious limit. This does not contradict the known fact that there 
are special nerve fibers for conduction of pain. Let P;, denote 
the sensitivity to a stimulus of modality «;, while Psa s¢B) ‘denotes 
the sensitivity to the stimulus of modality a, and intensity Bs 
Then we have 

Ps x58) © Psa, (3) 
and the sensitivity to pain caused by stimuli of modality a, is 
represented by 


PS asem = U Ps 4; ¢B) (4) 
B>B 
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where f, is the threshold of pain. The sensitivity to painful stim- 
uli of any modality is then represented by the subset of Py: 


Psa a LJ Ps sem)" (5) 


Psa; CPs 


The set which represents the property of conduction, Pc, con- 
sists of subsets P., which are conductions at different rates. 
Other subsets are conductions mediated chemically and conduc- 
tions mediated electrophysiologically. 

Of special interest is the set Py of all orderly movements. In 
this set are included such movements like those of cilia and flagella, 
as well as the movement involved in the protrusion of a pseudopod. 
The fact that different pseudopoda may succeed each other in a 
rather random fashion does not make the movements in each pseu- 
dopod disorderly or random. On the other hand, movements of trans- 
port due to either diffusion or to active transport, secretion or ab- 
sorption, are also su»sets of P,. Thus one important subset of 
Py is the set P, of all movements which may be called movements 
on the physically molar level. Another important set is the set of 
orderly movements P;,, on the molecular level. 

The set P, of properties of absorption includes as subsets the 
absorption Ge, of building substances for the organism and the 
absorption Bas of energy-yielding substances. The absorption P, f 
of energy is an element of P,. That absorption may either proceed 
through absorption’ Pa, of radiant energy directly, like in auto- 
trophic plants, or through absorption of energy-yielding substances, 


Pay We have 
OLS and Lies (6) 
Hence 
Py Pa, #9- (7) 


This gives us another example of (2). 

The set P,, of the properties of metabolism includes as subsets 
the anabolism or synthesis Py,,, catabolism Py, and storage Py, s 
Those contain other subsets. Thus, for example, synthesis of 
proteins is included in Pues} storage of fats and storage of sugars 
are both included in Raves , ete. 
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How are the different sets, P;, Py, etc., related in the primordial 
organism? We do not know the full answer to that question, but 
certain partial relations may be inferred from the logical analysis 
of what is usually meant by the concept of organism. An organism 
is a set S, of sets P, in which certain relations of ‘*immediate suc- 
cession’’ are established. Therefore those relations between the 
sets P, of S, can be best represented by a directed graph (cf. I). 

We define one biological property P, as immediately succeeding 
another, P,, if the manifestation of P, causes the manifestation of 
P, without involving the manifestation of any third property P, 
(1 4#i;1#k). As examples we may give excitation and conduction. 
However, to the extent that both P, and P, contain subsets P;, and 
Pg» a sequence of manifestations of different subsets may be in- 
volved in the sequence P,—» P,. Thus, for example, the excita- 
tion of a nerve is immediately followed by conduction of an im- 
pulse. But in different types of nerves both the excitation and the 
conduction may represent a sequence of somewhat different physi- 
cochemical processes. If P, succeeds P;, then we say that P, pre- 
ceeds P, . 

There will be a general agreement at least on the structure of 
some parts of the primordial graph. Thus we have in general 
P; — P,, and P, — P,. In some organisms we have only Puy: 
while in others both Puy and Pu, are present, but at least one of 
the two sets of properties is exhibited by a// organisms. As a re- 
sult of Py, as well as of Py, Wwe have in general ingestion, P, , or 
excretion of substances. Thus P, — P,. The process of diges- 
tion, which corresponds to a set P, of properties, occurs either as 
a result of Py, or of Py, or both. Hence we also have Py — P). 
Then again P, — Py. Metabolic processes being essential to the 
proper performance of sensitivity, conduction, and motion, we have 
also Py, — Ps; Py. — Pc; Py. — Py- The metabolic processes 
are also essential for the reproduction P, of the organism. Hence 
Pye — Pr, and by the same token, Py — Pp. 

There is even in the primordial organism also a set P,, of prop- 
erties of homeostasis. The connection between P, and the other 
P, is not directly evident from @ priori considerations. However, 
undoubtedly there is a connection Pye — Py, a8 well as a con- 
nection of P,, to nearly all other P,’s. 

We have already remarked that in a multicellular organism the 
different P,, (% = 1,2,8...) map on the P, of the primordial. We 
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have studied previously (I,V) different kinds of possible mappings. 
We shall consider here the most general conceivable mapping. Let 
S% be the set of sets P, which represents a multicellular organism, 
and S;—the set of sets P, which represents the primordial organism. 

Those sets are not yet spaces, since no topology has been in- 
troduced into them. By connecting the different points P, of Sp or 
P,., of S, by appropriate arrows, which indicate the relations of im- 
mediate succession, we obtain a one-dimensional continuous space, 
if we consider the arrows as continuous lines or as one-dimensional 
simplexes. But actually those arrows merely indicate the relation 
of immediate succession. As simplexes they have no physical 
meaning. It is therefore much more logical to construct spaces S, 
and S, by introducing into the sets Sj and S; a topology by an ap- 
propriate definition of neighborhoods. Considering the P,,, and P, 
as points in the topological spaces S$, and S,, respectively, this 
this can be done in the simplest manner by defining the neighbor- 
hoods in S, and S, in the following way: 

The neighborhood of a point in Sy (or S,) consists of the point 
itself and of all those points which it immediately precedes. 

Since the relation of immediate succession can be represented by 
an arrow of a graph, we see that to each S, and to Sp, there isa 
corresponding graph. Instead of saying that a point a immediately 
precedes 5, we may say that ‘‘a is connected to 6 by an arrow.”’ 
Then the above definition of neighborhoods can be stated thus: the 
neighborhood of a point in S, (or S,) consists of the point itself 
and of all points to which arrows originate from it. 

We shall subject the mapping of S, onto S, only to the condition 
of being continuous. If now the graph of the primordial organism is 
given, the possible spaces S, are obtained in the following manner: 

As we have remarked above (p. 73), the more complex an or- 
ganism the ‘‘narrower’’ the specialization of its cells; in other 
words, the smaller a subset P; that each cell type represents and 
therefore the larger the number n of subsets. Denote by P{” any 
set of n subsets of P,. Denote by Pra (r = 1,2,3,...,n) the n sub- 
sets of P{™. For example, if i =, and n = 2, and if P, had only 
three possible subsets, P;,, Ps., Psg; which denote respectively 
the sensitivities to light, sound, and touch, then Pee would be 
either (P,,, Psg), or (Psy Ps 3), OF (Psa P;,)- Actually, of course, 
each P; includes much more than two subsets. We now proceed as 
follows: 


78 N. RASHEVSKY 


(A). For each set P, in Sp, which is represented by a point on 
the corresponding graph, we choose n, subsets Pini . These sub- 
sets, for all possible values of 7, will be the points of the space 
Sy- Ifin S, the point P; is sonmected by an arrow to the point P, , 
then in S, at least one point Be "? is connected by an arrow to at 
least one aay: k, the directions “of the arrows being properly pre- 
served. If in S,, P,—» P,, then no P,,, is connected to any Py ge 
(a means ‘‘not a’’.) 


In this manner neighborhoods of S, are mapped on neighborhoods 
of Sp. The mapping is continuous. 

The proposed representation of S, and S, as topological spaces, 
in which neighborhoods are defined in the above manner, is more 
consistent than their representation by one-dimensional complexes 
or graphs. In the graphs used in I-VI the directed lines, or arrows, 
stand only to represent certain relations between the biological 
properties, namely, relations of immediate precedence and immedi- 
ate succession. As one-dimensional complexes those lines have 
no physical meaning. To the extent, however, that to each Sj), as 
well as to the S,, there corresponds a graph, we may conveniently 
use the terminology of graphs in our discussion. But actually 
while a graph, considered as a one-dimensional complex, is a 
polyhedron in a metric space, S, and S, are topological spaces, 
which are not metrized up to this point. This point of view is in 
line with Alexandroff’s and Hopf’s (19385) attitude of not separating 
the combinatorial and set-theoretical methods in topology. 

The above procedure of constructing S, from S, does not prevent 
tie spate that see than one pin will i Connected to a 

Pe , or that one pir is connected to several P, “. The total 
number of possibilities for prescribed values of n; (¢ = 8,C,M, etc) 
is very large. To each possible set of connections permitted by 
(A) there corresponds a different multicellular organism. Moreover, 
it must be apmomparsq that the set P, can be divided in many dif- 
ferent ways into oS subsets. Thus we may have two sets S, 
which are boikectinie but in which two topologically identical 
points Prat Pil represent. different subsets of P;. Such two sets repre- 
sent again ‘different multicellular organisms. 

It is an interesting combinatorial problem to determine the total 
number N of possible different sets S, which correspond to a given 
Sp. This number will give us the total number of possible or- 
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ganisms. If any of the P.’s are infinite sets, then of course N will 
be infinite. As we femtinrkedl actually for any given organism, n; 
are always finite. But inasmuch as theoretically some, if not at 
of the a s may be infinite and even not enumerably infinite, we 
have here a possibly interesting problem in set theory. Namely, 
for given cardinal numbers of the P,’s what is the cardinal number 
of the set of possible sets S,? 

The organism may be eee as the more complex, oe more 
connecting arrows there are in S, between the different Ps "i and 
ee : the maximum degree of aamiees is obtained when each 
point oo "* is connected by an arrow to each oe: We shall call 
such a S, the maximal S, 

The ate (A) for constructing S, from S, is very general and con- 
tains only the restriction of continuity. We shall, however, add to 
it another restriction. A space S,, constructed according to (A), is 
not necessarily connected. To see this we may consider the par- 
ticular case when all n,’s are equal (n; =n). We may then connect 
each point Pe of S, to just one corresponding point PY in such a 
manner that we have n disjoined graphs, each of which is homeo- 
morph with S,. It is to be noted that with the above definition of 
neighborhoods in S, a connected graph implies a connected space 
S,, in the more general topological meaning of the word ‘*connec- 
tion.’’ Biologically this would mean that a multicellular organism 
simply consists of n separate organisms, which in no way interact 
with each other. To avoid this absurdity we add the requirement: 


(B). S, is connected. 


Requirement (B) implies some geometrical and combinatorial re- 
lations in S,, relations which are in general by no means simple, 
but which may have definite biological predictive value. We shall 
not discuss here this question in detail, but merely illustrate its 
nature on an example. 

Let again n, =n, and consider the graph of So. In the absence 
of requirement (B) this graph may consist, as we have seen, of n 
disjoined identical graphs. 

Because of n, =n, we can conveniently nuniber all the Pea and 


P2> in such a manner that g=7, in other words, so that Ba 
c 


80 N. RASHEVSKY 


To make the whole graph connected, we must have at least n — 1 
lines joining different non-corresponding points of the n identical 
graphs. Let such a cross-connection be made between the points 


(n) (n) : (n) (n) iS 
Pre, and Perce yf Since we already have Pins —_ Ei tes » there 


fore the point rad will have two predecessors, namely, Pie and 


ey Each of the last two points will have at least one im- 
mediate predecessor. 

We call a point a a predecessor of b, if a precedes 5. However, 
a may precede 6 either immediately by being connected to 6 by an 
arrow a — 6b, or not immediately by being merely connected to 6 by 
a sequence of arrows, or by a directed way (III, p. 114). Accord- 
ingly we shall speak of immediate and non-immediate predecessors. 
Hence, requirement (B) assures us that of the n points which corre- 
spond to P,, n —1 points will have more predecessors in S, than 
P, has in Sp. In other words, while a certain point P, in Sp may be 
reached from only one point P,, in S, we shall have at least n — 1 
points corresponding to P,, each of which will have at least two 
immediate or non-immediate predecessors. The maximal S, is of 
course always connected. 

If P; is a predecessor of P,, then anything that affects P. is 
going to affect P,. Hence the biological meaning of the above is 
that in a multicellular organism, due to requirement (B) the number 
of biological properties which affect a given property is greater 
than in a unicellular organism. 

As in I, we prove in quite a similar manner the following: 

Theorem I. The rules (A), when applied to two different primordial 
graphs S, result always in different graphs S,. 

On the other hand, the theorems in IV require some modifica- 
tion. The general situation is now much more complicated. For 
the case of a maximal S, we have the following: 

Theorem Il. If in the graph of S, there is a way (III, p. 114) 
from P; to P, , which goes through one or more intermediate points 
P,, then in the maximal S, the number of disjoined ways from any 
point Pig to any pyre is equal ton, .., where n_,. is the smallest 
of all n,’s, which belong to any of the intermediate points P.. 

Proof. If n,,, is the smallest of all n,’s, then we can pick out 
Nin POints amongst every of the n,’s, and construct Nmin disjoined 


(n;) (n,) woes : 
ways from Prot, to P, a. - More than n, disjoined ways are im- 
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possible, because some of them would have common points amongst 
the n_,,, points of Bote, 

How about unicellular organisms? They exhibit a great variety 
of properties P;, and they also map on each other and on the pri- 
mordial. Some microorganisms, as we remarked above (p. 76), do 
not have any Py, (non-motile microorganisms), but all of them have 
Puyo: Similarly different microorganisms have different Psy» but at 
least some elements of P, occur in all of them. Generalizing by 
induction, it is natural to postulate: 


(C). The spaces S,, which represent all possible unicellular or- 
ganisms are obtained from S, by assigning to one cell any possible 
space S,, with the exception of such elements Px or combina- 
tions thereof as would be incompatible with the size of a cell or 
with physical laws. 


The last restriction is essential for example in the following 
case: Consider the property P, of homeostasis. Amongst its ele- 
ments the set P, contains the subset P, of temperature regulation. 
While some microorganisms show other homeostatic properties P,, , 
for example osmoregulation (Prosser, 1950), no unicellular organism 
can show temperature regulation, because, as simple calculations 
show, the maintenance of any appreciable difference in temperature 
between a single cell and its environment would require such rates 
of heat producing reactions as are incompatible with known physi- 
cochemical laws, because of the tremendous specific surface of a 
single cell. 

In addition to (C), we postulate: 


(D). All possible multicellular organisms are represented by 
S,), and conversely, any organism which corresponds to an S, is 
possible. 


A number of interesting conclusions follow almost immediately 
from (A)-(D), without a more detailed specification of the structure 
of the primordial space S,.- 

From (A) it follows that if we have an S,, in which P,, — P,, 
then in the primordial we have P, —» P,. Hence, from (D) it follows 
that there exist organisms, in which P;,, — P,,, where p and vare 
in.general different from « and 8. From (C) it follows that there 
exist unicellular organisms, in which P,,, — P,,, unless the exist- 
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ence of either P;,, or P,,, or both, is incompatible with the size of 
a single cell. 

Thus from observations of some relations within an organism, we 
may infer the existence of other organisms in which different but 
corresponding relations exist. 

Let us consider the sequence P; —> Pp — Py. As we have 
seen, an important subset of Py is Py,, the subset of molar orderly 
movements. In a unicellular organism this may be represented by 
amoeboid movements or ciliary or flagellary movements, which we 
shall denote correspondingly by Pua? Pus Puy: We have 


Pye CPt PRC Py Gamer eth Re (8) 


In a higher animal different cells specialize in the three different 
subsets; the leucocytes, for example, in Py’ cells of respiratory 
epithelium in Pui spermatozoids in Py, 

A subset of Puy is represented by contractile muscular tissues. 

A movement of a unicellular organism, for example an amoeboid 
movement, plays several roles in the life of the organism. 

(a) It propels the organism as a whole in space, towards food or 
away from the enemy; 

(b) it orients the organism differently with respect to different 
stimuli, thus 

(c) changing possibly the effectiveness of the stimulus; 

(d) it helps ingestion of food into the digestive vacuole, and 

(e) it may cause an eventual movement of food in the latter. It 
also 

(f) in general affects the transport of substances inside the cell. 

The sets of movements which produce correspondingly the ef- 
fects (a), (b), (c), (d), (e), and (f) are included in Puy: In line with 
(A), (B), and (D), different cell types specialize in some higher 
organisms in each of those subsets of movements, which we shall 


denote correspondingly by Pca); Pucsy Pycgs etc. We have 


Puca) C Py (9) 


and similarly for Py,,, and the others. 

The skeletal muscles are general representatives of Pyca)s they 
also represent the subset Parcs)’ The ciliary muscle of the eye or 
the muscles of the tympanic membrane have properties which be- 
long to Pyc,+ 
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Instead of the general sensitivity P, we also have in a higher 
animal several subsets of P;, some of which we discussed above. 
The same holds for P, 

Let Ps, denote the set of visual sensitivities of all kind. We 
have 


Bake a te (10) 


In a primordial we have P,; —> P. —> Py. Hence it follows from 
(A) that in some higher organisms we must have: 


Ps, > Pe > MCd)* (11) 


In words: There exists higher organisms, in which gastrointestinal 
movements are affected by visual stimuli. A well-known example 
is the vomiting at the sight of some unpleasant things. By the 
same argument, we have more generally 


a me ae (12) 


which states that in some animals different sensory stimuli affect 
the gastrointestinal motility. 

As we have seen, the more complex an animal, the more subsets 
P;,, are connected to a given subset P,,, if P,—» P,. Hence, the 
higher the animal, the greater the variety of sensory stimuli, which 
affect the gastrointestinal motility. Man is the highest of the 
presently existing animals, though there is no reason to assume 
that it is represented by a maximal S,. We should, however, expect 
the human gastrointestinal tract to be more sensitive to a great 
variety of sensory stimuli, than in any of the other animals. 

To the extent that our psychological experiences are manifesta- 
tions of a complex interplay of sensory activities, we may state 
that it follows from (A), (B), and (D) that in man, as well as in 
some animals, psychic disturbances produce gastrointestinal trou- 
bles. This is a well-known fact. This does not follow, however, 
and cannot follow from any metric biophysical theory, unless spe- 
cifically assumed. It does follow, however, from the general 
principle of biotopological mapping, as formulated here. It cannot 
be deduced from considerations of natural selection, since its 
adaptive value is, if anything, negative, at least for the usual 
occupations of man. 

By the same argument we see that there exist animals, in which 
sensory stimuli and specifically in man, psychic disturbances, 
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affect the cardiovascular phenomena. We merely substitute in (12) 
Pycsy for Py cay In this connection it is interesting to note that 
‘“‘stimulation of almost any afferent nerve of the body can affect 
the heart rate’’ (Carlson and Johnson, 1953, p. 161). 

From P, — P, — Py it follows, according to (A), that 


Pou aot Care Puy» (13) 


where P,., is any subset included in P, - 

Hence we may substitute for Py, any subset contained in Pus: 
But Py, represents the different molecular orderly movements, such 
as, secretion. In particular the set Py, of ‘secretions of different 
digestive enzymes’’ is included in Pup and hence in Py. Therefore 


we have 
Ps.— Pog— Py, - (14) 


In words: In some animals sensory stimuli affect the secretion of 
digestive enzymes. In men we should expect psychological states 
to affect such a secretion. 

A biological property of an organism is affected whenever any 
predecessor is affected. If a given P., has, say, n predecessors 
and if we assume equal probability for any of them being disturbed, 
then the probability that an observed disturbance of the property 
P,;,, is actually due to a disturbance of a predecessor, and not of 
P,,, itself, is (1 -1)/n. To know the total number of predecessors 
for a given P,;, we must know the structure of Sp, which we do not 
know. But even from a simple general consideration as that made on 
page 76 we see that for P,,, and therefore for any Pucay OF Pys 
n> 4, at the very least. Hence the probability of a given gastro- 
intestinal disturbance being due to disturbances of other parts of 
the organism is greater than 0.75. To the extent that the psycho- 
logical disturbances involve P,.’s and Pcog’s, we can infer that 
the percentage of all observed gastrointestinal disturbances, which 
is due to psychic disturbances, is greater than 60%. It is said to 
be actually over 80%. Consider now in Sp the relation of Py to 
some following P,’s. Molar movements are followed by feeding re- 
actions P;. So are the molecular orderly movements, because 
diffusion gradients and active transport bring food from the sur- 
roundings towards the cell in non-motile cells. Hence 


Py — Pr (15) 
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in the primordial. Denote by Pucsec) the subset of Py which con- 
tains all the secretory properties. Since Pucsec) © Py» therefore 
from (15) and from (A) and (B) it follows: 


E (S,): Euitabaxeeoet resol Stray: (16) 


(sec) 


In words: There are organisms in which secretory phenomena 
produce feeding reactions. An example of this is offered by Urechis 
which secretes mucous bags that filter out food particles and are 
eventually swallowed (Prosser, Joc. cit. p. 145). Another example 
is the secreted spider-web which catches the prey. 

In some unicellulars the movements P,,, affect the transport of 
substances both within the cell (by stirring ‘‘the protoplasm’’) and 
outside of it. (Ingestive movements of cilia.) Denoting phenomena 
of transport in general by P,, we thus have: 


E (S94)? Pyy © Soy-Pr © Soy-Py; — Pr (17) 
Hence it follows from (A) and (C): 
Beets PtS 5. 96P acct: Ps (18) 


where > is the sign of implication. Hence, in the primordial we 
have P,, —> P,. But then it follows from (A): 
E (So): Pysa © So-Prg © So-P. 


win Pre: (19) 


Put for P,,,,, the movements of skeletal muscles, which are pro- 
pelling the organism in space. Put for Pp the set of circulatory 
phenomena. Then (19) states that in some higher organisms such 
movements of skeletal muscles affect the circulation. Example: 
The venous circulation in man is largely due to the ‘‘kneading”’ 
action of skeletal muscles. Effects of body musculature on circu- 
lation are also known in a number of lower animals (Prosser, loc. 
Cit.). 

The general principles (A)-(D) apply to both animals and plants. 
The latter are perhaps characterized by the fact that ee = 0; the 
subset of all amoeboidal movements of the cell is empty. We have 
Pu, © Pus: But Py, # 0, since we have the flow of water in plants. 
The set of flagellate movements is not empty in lower unicellular 
plants. Similarly the subset of P, which contains the conduction 
due to electrochemical phenomena (nerve conduction) is empty in 
plants. However, from P; — P, — Py in the primordial it follows: 
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This applies to plants as well as to animals. If P;,, denotes any 
sensitivity to external stimuli, P;-g—chemical mediation, and Py. 
denotes either the movement of water through tracheae or the move- 
ment. of food-containing solvents through sieve cells, then (20) 
states that adequate stimulation of plants, say by light, affects 
both the flow of sap and the translocation of metabolites. 

The subset of all secretions, which is included in P,,,, includes 
such subsets as secretion of wastes, secretion of poisons, secre- 
tion of substances useful to other organisms, and secretion of sub- 
stances lethal or otherwise harmful to other organisms. From (C) 
it follows that there must exist microorganisms which secrete sub- 
stances that are either lethal or otherwise harmful to some other 
organisms. The toxins of some pathogenic bacteria is one ex- 
ample. The set of substances harmful to some other organisms in- 
cludes the subset of substances harmful to some microorganisms. 
The antibiotics, like penicillin, are an example. 

By the same token there must exist microorganisms which pro- 
duce substances that are useful to other organisms, which is a well 
known fact. 

The finding of different possible subsets of a set P,; may be a 
matter of logic. A set P, consists of all conceivable subsets in- 
cluded in it, as long as the existence of such subsets is com- 
patible with physical laws. Sometimes the possibility of a subset 
P,,, is shown by experiment. 

In many relatively higher animals there are cells which upon 
stimulation with appropriate substances, the antigens, produce 
special proteins, the antibodies. The production of antibodies is 
thus a subset of the set P,,, of metabolic phenomena. The sen- 
sitivity to antigen is a subset of the set P;. We shall designate 
the former by Py.,q), the other by P5,,). From (C) it follows that 
there must exist unicellular organisms which produce antibodies 
when stimulated by an appropriate antigen. So far no such micro- 
organism has been definitely found, though there is a report of 
antibody production by yeast (Fox and Plaisted, 1953; Cushing and 
Campbell, 1957). 

This is an example of prediction made from the general principles 
proposed here. 
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Another conclusion is also interesting. We have Ps) © P, and 
Pyeca) © Pye: We also have in higher organisms: 


ey es ees eh eee de ore (21) 
where P,,—>... Py denotes the possible intermediate sequence 


of points of S$), which as yet are unknown. From (A) and (21) it 
follows that in a primordial 


eal eS Ne Sw ae (22) 


and then again from (A) that there exist multicellular organisms in 
which 
9S ima ene eee ae (23) 


m@ Me(a)°* 


Here Po, stands for sensitivity to any stimulus, such as light, 
heat, etc. Expression (23) states that in some cases antibodies 
may be formed in multicellular organisms by ‘‘physical stimuli.” 
This reminds us of the “‘physical allergy’’ phenomena (Boyd, 1956), 
which seem, at least in some cases, to be connected with actual 
formation of ‘‘something’’ that can be transferred passively by the 
serum (Sherman and Seebohm, 1950). Alternative explanations of 
the observed phenomena seem, however, possible (Boyd, loc. cit.). 
Thus the above conclusion cannot be regarded as verified. Neither 
does it seem to be yet disproved. 

By the same line of reasoning we should expect that some plants 
produce antibodies. Indications of the existence of such a phe- 
nomenon have been found (Wallace, 1950). Again, however, a 
different interpretation of Wallace’s observation has been suggested 
(Luria, 1950). The evidence against production of antibodies by 
plants is considered by some to be conclusive (Cushing and Camp- 
bell, loc. cit.). However, in view of the fact that still very few 
plants have been investigated, the search would be worth while 
to be continued. Should, however, such a search fail definitely, 
this would lead to definite conclusions which are again verifiable 
by experiments. As has been remarked, in plants the set Py,, of 
amoeboid movements is empty. There are other subsets P;,, which 
are empty in plants. Plants are characterized by the set P, of 
subsets Pj, which are empty in them. If plants do not produce 
antibodies, then Py,,,) is empty in plants. Either then Py,(,) is an 
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element of P,, and in a sense enters into the definition of plants, 
or there must be some physicochemical relation between some of 
the P;, and Py,,,), such that the absence of Py, implies the ab- 
sence of Py,,,)- The first possibility does not seem very likely. 
The second possibility suggests immediately a possible physico- 
chemical connection between Pye and Py,¢q)- There may be a 
common biophysical or biochemical factor involved in certain type of 
movements of the cell and production of antibodies. It may, per- 
haps, be significant that antibodies are apparently produced by 
reticuloendothelial cells which resemble more than others some 
amoeboid forms. If there is a relation between Py, and Pyecay> 
then, however, we should look for antibody production by unicellu- 
lars not to yeasts, but rather to amoebae. 

The above is meant to illustrate that even in its immediate and 
simplest consequences, the postulates (A)-(D) do have a predictive 
and heuristic value. More of this should be found by further elab- 
orating the set-theoretical properties of S$, and S,. 

We have considered two subsets of ee the subset of what may 
be called nervous conduction and the subset of humoral chemical 
conduction. We shall denote the former by Po,- The set Po, in- 


cludes two important subsets, Pet and P-,,, such that 


Poig C Pe, and Po, © Pe: (24) 


The first one, Pei? we shall call, for lack of a better term, non- 
hysteresis conduction; the second, Poi,» hysteresis conduction. 
The first is characterized by independence of a given conduction 
phenomenon on the past history. The other is characterized by the 
fact that the character of the conduction depends on past con- 
ductions. No more detailed specifications are necessary here. In 
its simplest form a hysteresis conduction may manifest itself 
merely by the fact that a given response to a given stimulus is 
facilitated by repetition. In its most complex manifestations hys- 
teresis conduction may result in most complicated phenomena of 
learning in higher animals. Whatever physiological or biological 
theory of learning we accept, whether we consider it based on the 
existence of self-circuited neurons (Rashevsky, 1948) or on synaptic 
changes of a physical or chemical nature (Shimbel, 1950), or any 
other different conceivable mechanism, in all of them we deal with 
hysteresis conduction in the central nervous system. Conduction 
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here is understood not only along a nerve-fiber, but including syn- 
aptic transmission as well. 

From (C) it follows that there must exist unicellular organisms 
which exhibit some form of hysteresis conduction. The turning 
around of a paramecium in a capillary tube with fewer turns after 
several trials is an example of this (Prosser, loc. cit., p. 842). 
Whether we call such phenomena elementary acts of learning or not 
is immaterial. Relationally they are isomorphic to the more com- 
plex phenomena of learning in higher animals. 

So far as relations between the different P,’s and P,,,’s are con- 
cerned, we have considered only the relation of immediate preced- 
ence FR, and immediate succession FR, its converse. If we represent 
S, and S, by directed graphs then we have here a graphical repre- 
sentation of the class of relations R. We can apply to this some 
standard expressions of the theory of relations, but we do obtain 
only relatively trivial results. We give here just two illustrations. 

Thus, if we use the rather convenient notations of J. Riguet 
(1948) and denote by pr, the argument of the relation 2, we have 
the formula: 


pr, SE Cpr, Rk, (25) 


where SR denotes the product of the relations S and Fr. 
In particular we have 
pr, R* Cpr, R. (26) 


If R is the relation between P,,, and P,g when P,,, immediately 
precedes P, g, then R? is the relation between two points such that 
there is a third one, immediately preceding the second one and 
immediately succeeding the first. We shall say that the second 
point is a second-order successor of the first one. In that case 
(26) states that the number of points P, which have immediate suc- 
cessors is less or equal to the number of points which have second- 
order successors. 

Consider the property: ‘‘Cells are sensitive to different forms of 
light’’ as a relation between the cell c and light 7. In symbols cl. 

Let in the relation zy, R(x) denote the set of all y’s, such that 
aRy. Let 


R(xX)=[(_] R@;  RIXl = [| Re). (27) 


x EX xEX 
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Then (Riguet, Joc. cit.) we have: 
XCXe DRUG WER Xone (28) 
RO OXe sD Rid, | CREM (29) 


R(c) is the set of all forms of light to which a given cell is 
sensitive. If X is a given set of cells, then R(X) is the set of all 
forms of light to which we find sensitivity amongst cells of the 
set X; while R[X] is the set of all forms of light, such that every 
cell of X is sensitive to all of them. If X, is the set ofall cells 
of an organism and X, a set of cells of a part of it, then expres- 
sions (28) state that the number of different forms of light to which 
we find sensitivity amongst the cells of a part of the organism is 
less or equal to the number of all forms of light to which we find 
sensitivity amongst cells of the whole organism. Expression (29) 
states that the number of different forms of light such that every 
cell of a given part of an organism is sensitive to it is greater or 
equal to the number of different forms of light such that every cell 
of the organism is sensitive to all of them. 

By applying expressions (28) and (29) to the inverse relation R 
(different forms of light affect different cells), we find even more 
trivial statements: 

The number of cells which are sensitive to any wave length of 
light is greater or equal to the number of cells which are sensitive 
to only some selected spectral range. And: The number of cells 
which are sensitive to all wave lengths within a given spectral 
range is less than or equal to the number of cells which are sensi- 
tive to all wave lengths within a wider spectral range. 

The above statements while true have no predictive or heuristic 
value. It is, however, rather likely that by studying in a similar 
manner more complicated relations between pairs of points in So 
and S, we shall arrive at much more valuable conclusions. 

A systematic application of theory of relations to biology and an 
attempt to build an axiomatic biology on this basis has been made 
in the noteworthy research of J. H. Woodger (1937). Though 
Woodger’s approach differs seemingly very radically from ours, the 
possibility cannot be denied that the further development of both 
approaches may establish a number of points of contact between 
them. In any case, to Woodger belongs the credit of having clearly 
emphasized the importance of relational aspects in biology and to 
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have made the first systematic mathematical study of those rela- 
tional aspects. 

The principal difficulty found in the present approach is that the 
graphs of Sp, or S, do not represent merely such simple relations 
as immediate succession. Two biological properties of which one 
immediately succeeds the other are usually characterized by other 
relations which are more inclusive than the relation R of immediate 
succession. 

No elaborate scientific system can be developed on the basis of 
one or two general principles only. Other principles or postulates 
will have to be added in the future to those presented here. Some 
of them are likely to be of a restricting nature, reducing somewhat 
the generality of (A) or (C), just as (B) does restrict (A). As we 
have already suggested in I, the principle of maximum simplicity, 
introduced by us (Rashevsky, 1948), or as David Cohn (1954, 1955) 
aptly calls it, the principle of optimal- design, may be used as 
another general mathematical principle, together with the principle 
of biotopological mapping. The principle of optimal design may 
possibly eliminate some of the possibilities presented by (A). 

The general formulation presented here, unlike the ones suggested 
in I-VI, does not give any indications as to arrangements of vari- 
ous differentiated cells into organs. It seems that an additional 
rule of a relational character must be added to (A)-(D). 

It may be asked as to whether we should not expect on the basis 
of the foregoing to have organisms which are sensitive to X-rays 
or y-rays, or which can perform other tasks than those so far ob- 
served. The answer is in the affirmative, since, for example, 
sensitivity to y-rays is included in sensitivity to radiation. Does 
it then follow that we must predict organisms to develop in the 
future, which will have such properties? Not necessarily. 

Through science and technological invention man has found 
means to detect y-rays and to perform numerous tasks which no 
organism can do directly. All such performances we do not con- 
sider as part of the biological properties of the human organism. 
But all such performances are definitely the result of biological 
manifestations of the human organisms, in particular of the brain. 
And relationally it is quite consistent to extend the notion of or- 
ganism to include the results of the direct biological properties. 
Relationally there is no difference between the Urechis secreting 
a mucous bag which catches food and then swallowing it, or the 
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spider secreting its web, and a human-being manufacturing as a re- 
sult of his brain-work a machine which processes food and makes 
it possible for it to be swallowed. 

Such an extension of the relational principles will necessitate 
a more detailed study of relations between an organism and its 
surroundings, which thus far have been considered only sketchily. 
The principle of optimal design may then lead to the conclusion 
that organisms which perceive X-rays or y-rays directly may not 
develop, because it is simpler to develop an appropriate brain 
which can invent technological devices for the performance of such 
tasks than to develop an organism which performs them directly. 

An extension of the relational principles, such as mentioned 
above, will also include many aspects of social relations, both in 
man and animals. 


The author is indebted to Mr. Robert Rosen for checking the 
manuscript. 

This work was aided in part by a grant from the Dr. Wallace C. 
and Clara A. Abbott Memorial Fund of The University of Chicago 
and in part by United States Health Service Grant RG~5181. 
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ANNOUNCEMENT BY PERSPECTIVES IN 
BIOLOGY AND MEDICINE 


PERSPECTIVES IN BIOLOGY AND MEDICINE, a new quarterly 
journal dedicated to a multidisciplined approach to the problems of 
biology and medicine, is announced by the Division of Biological 
Sciences and the Press of The University of Chicago. 

PERSPECTIVES will present new hypotheses and concepts 
representing informed thinking; interpretive, selective essays which 
take stock of the implications of recent and current research and 
indicate strategy for the future; and, because perspective and in- 
spiration can also be gained by looking to the past, reviews by 
great men in science and medicine of their formative years and 
their philosophy of research. 


95 


